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EDITOR’S PREFACE 

The editors believe that the reader who has finished the study of this 
book will see the full justification for including it in a series of volumes 
dealing with aeronautical subjects. 

However, the editor’s preface usually is addressed to the reader who 
starts with the reading of the volume, and therefore a few words on our 
reasons for including Professor Michal’s book on matrices and tensors 
in the GALCIT series seem to be appropriate. 

Since the beginnings of the modern age of the aeronautical sciences 
a close cooperation has existed between applied mathematics and 
aeronautics. Engineers at large have always appreciated the help of 
applied mathematics in furnishing them practical methods for numerical 
and graphical solutions of algebraic and differential equations. How- 
ever, aeronautical and also electrical engineers are faced with problems 
reaching much further into several domains of modern mathematics. 
As a matter of fact, these branches of engineering science have often 
exerted an inspiring influence on the development of novel methods in 
applied mathematics. 

One branch of applied mathematics which fits especially the needs 
of the scientific aeronautical engineer is the matrix and tensor calculus. 
The matrix operations represent a powerful method for the solution of 
problems dealing with mechanical systems with a certain number of 
degrees of freedom. The tensor calculus gives admirable insight into 
complex problems of the mechanics of continuous media, the mechanics 
of fluids, and elastic and plastic media. 

Professor Michal’s course on the subject given in the frame of the 
war-training program on engineering science and management has 
found a surprisingly favorable response among engineers of the aero- 
nautical industry in the Southern Californian region. The editors be- 
lieve that the engineers throughout the country will welcome a book 
which skillfully unites exact and clear presentation of mathematical 
statements with fitness for immediate practical applications. 

Theodore von KArmAn 
Clark B. Millikan 
vii 




PREFACE 


This volume is based on a series of lectures on matrix calculus and 
tensor calculus, and their applications, given under the sponsorship 
of the Engineering, Science, and Management War Training (ESMWT) 
program, from August 1942 to March 1943. The group taking the 
course included a considerable number of outstanding research en- 
gineers and directors of engineering research and development. I am 
very grateful to these men who welcomed me and by their interest 
in my lectures encouraged me. 

The purpose of this book is to give the reader a working knowledge 
of the fundamentals of matrix calculus and tensor calculus, which he 
may apply to his own field. Mathematicians, physicists, meteorologists, 
and electrical engineers, as well as mechanical and aeronautical en- 
gineers, will discover principles applicable to their respective fields. 
The last group, for instance, will find material on vibrations, aircraft 
flutter, elasticity, hydrodynamics, and fluid mechanics. 

The book is divided into two independent parts, the first dealing 
with the matrix calculus and its applications, the second with the 
tensor calculus and its applications. The minimum of mathematical 
concepts is presented in the introduction to each part, the more ad- 
vanced mathematical ideas being developed as they are needed in 
connection with the applications in the later chapters. 

The two-part division of the book is primarily due to the fact that 
matrix and tensor calculus are essentially two distinct mathematical 
studies. The matrix calculus is a purely analytic and algebraic sub- 
ject, whereas the tensor calculus is geometric, being connected with 
transformations of coordinates and other geometric concepts. A care- 
ful reading of the first chapter in each part of the book will clarify 
the meaning of the wprd “tensor,” which is occasionally misused in 
modern scientific and engineering literature. 

I wish to acknowledge wuth gratitude the kind cooperation of the 
Douglas Aircraft Company in making available some of its work in 
connection with the last part of Chapter 7 on aircraft flutter. It is a 
pleasure to thank several of my students, especially Dr. J. E. Lipp 
and Messrs. C. II. Putt and Paul Licber of the Douglas Aircraft 
Company, for making available the material worked out by Mr. Lieber 
and his research group. I am also very glad to thank the members of 
my seminar on applied mathematics at the California Institute for 
their helpful suggestions. I wish to make special mention of Dr. C. C- 
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Lin, who not only took an active part in the seminar but who also 
kindly consented to have his unpublished researches on some dramatic 
applications of the tensor calculus to boundary-layer theory in aero- 
nautics incorporated in Chapter 18. This furnishes an application of 
the Riemannian tensor calculus described in Chapter 17. I should 
like also to thank Dr. W. Z. Cliien for his timely help. 

I gratefully acknowledge the suggestions of my colleague Professor 
Clark B. Millikan concerning ways of making the book more useful 
to aeronautical engineers. 

Above all, I am indebted to my distinguished colleague and friend. 
Professor Theodore von Karman, director of the Guggenheim Graduate 
School of Aeronautics at the California Institute, for honoring me by 
an invitation to put my lecture notes in book form for publication in 
the GALCIT series. I have also the delightful privilege of expressing 
my indebtedness to Dr. Karman for his inspiring conversations and 
wise counsel on applied mathematics in general and this volume in 
particular, and for encouraging me to make contacts with the aircraft 
industry on an advanced mathematical level. 

I regret that, in order not to delay unduly the publication of this 
book, I am unable to include some of my more recent unpublished 
researches on the applications of the tensor calculus of curved infinite 
dimensional spaces to the vibrations of elastic beams and other elastic 
media. 


California Institute op TECHNOLoor 
October, 1946 


Aristotle D. Michal 
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PART /. MATRIX CALCULUS 
AND ITS APPLICATIONS 

CHAPTER 1 

ALGEBRAIC PRELIMINARIES 


Introduction. 

Although matrices have been investigated by mathematicians for al- 
most a century, their thoroughgoing application to physics/ 1 engineer- 
ing, and other subjects^ — such as cryptography, psychology, and 
educational and other statistical measurements — has taken place only 
since 1925. In particular, the use of matrices in aeronautical engi- 
neering in connection with small oscillations, aircraft flutter, and elastic 
deformations did not receive much attention before 1935. It is inter- 
esting to note that the only book on matrices with systematic chapters 
on the differential and integral calculus of matrices was written by 
three aeronautical engineers. J 


Definitions and Notations. 

A table of mn numbers, called elements, arranged in a rectangular 
array of m rows and n columns is called a matrix ^ with m rows and n 
columns. If a} is the element in the 2 th row and jth. column, then the 
matrix can be written down in the following pictorial form with the 
conventional double bar on each side. 


« 1 , 02 , ' • • , 
2 2 

Oi, O2, * ’ * , 


On 

On 


a 


m 

n 


In the expression a} the index i is called a superscript and the index j a 
subscript. It is to be emphasized that the superscript i in aj is not the 
2 th power of a variable aj. 

If the number m of rows is equal to the number n of columns, then 


t Superior numbers refer to the notes at the end of the book, 
t Frazer, Duncan, and Collar, Elementary Matrices and Some Applications to 
Dynamics and Differential Equations^ Cambridge University Press, 1938. 

1 
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the matrix is called a square matrix.^ The number of rows, or equiva- 
lently the number of columns, will be called the order of the square 
matrix. Besides square matrices, two other general types of matrices 
occur frequently. One is the row matrix 

II ai, 02, • • an II ; 

the other is the column matrix 

a2 


II II 

It is to be observed that the superscript 1 in the elements of the row 
matrix was omitted. Similarly the subscript 1 in the elements of the 
column matrix was also omitted. All this is done in the interest of 
brevity; the index notation is unnecessary when the index, whether a 
subscript or superscript, cannot have at least two values. 

It is often very convenient to have a more compact notation for 
matrices than the one just given. This compact notation is as follows: 
if is the element of a matrix in the ith row and jth column we can 
write simply 

II 4 II 

instead of stringing out all the mn elements of the matrix. In par- 
ticular, a row matrix with element a* in the fcth column will be written 

11 11 , 

and a column matrix with element aJ^ in the A;th row will be written 

II o* II • 


Elementary Operations on Matrices. 

Before we can use matrices effectively we must define the addition of 
matrices and the multiplication of matrices. The definitions are those 
that have been found most useful in the general theory and in the 
applications. 

Let A and B be matrices of the same type, i.e., matrices with the same 
number m of rows and the same number n of columns. Let 

A = II a} II . -B - II 4 II . 

Then by the sum A + B oi the matrices A and B we shall mean the 


t It will occasionally be convenient to write aij for the element in the 'ith row 
andyth column of a square matrix. See Chapter 5 and the following chapters. 
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uniquely obtainable matrix 

c^W 4 W, 

where 

c) = a)-\- bj (i ^ 1, 2, • • • , m; i = 1, 2, • • • , w). 

In other words, to add two matrices of the same type^ calculate the 
matrix whose elements are precisely the numerical sum of the cor- 
responding elements of the two given matrices. The addition of two 
matrices of different type has no meaning for us. 

To complete the preliminary definitions we must make clear what we 
mean when we say that ituo mairzces are cgwaZ. Two matrices .4 = || a} || 
and jB = II fey II of the same type are equal, written as A = i5, if and 
only if the numerical equalities a) = fey hold for each i and j. 


Then 


Exercise 



1, 

-1, 

V2, 

5 

A = 

0, 

0, 

3, 

-2 


1.1, 

2, 

-4, 

1 


0, 

0, 

-V2, 

1 

B = 

0, 

0, 

-1, 

3 


1, 

0, 

2, 

-4 


1, 

-1, 

0, 

6 

A+B = 

0 , 

0, 

2, 

1 


2.1, 

2, 

-2, 

-3 


The following results embodied in a theorem show that matric 
addition has some of the properties of numerical addition. 


Theorem. If A and B are any two matrices of the same type^ then 

A + B = B + A. 

If C is any third matrix of the same type as A and then 

(A + B) C = A + (B + C). 

Before we proceed with the definition of multiplication of matrices, 
a word or two must be said about two very important special square 
matrices. One is the zero matrix, i.e., a square matrix all of whose 
elements are zero, 

0, 0, --^O 

0 , 0 , - - • , 0 


0 , 0 , • • * , 0 
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We can denote the zero matrix by the capital letter O. Occasionally 
we shall use the terminology zero matrix for a non-square matrix with 
zero elements. 

The other is the unit matrix , i.e., a matrix 

/ - II II . 

where 

5* « 1 if 

= 0 if i ^ j. 

In the more explicit notation 

1 , 0 , 0 ,* • 0 
Oj > 0 

0, 0, 1, 0,* • - , 0 

• 


II 0,0, 1 II 

One of the most useful and simplifying conventions in all mathe- 
matics is the summation convention: the repetition of an index once as a 
subscript and once as a superscript will indicate a summation over the 
total range of that index. For example, if the range of the indices is 
1 to 5, then 

5 

af}^ means '^aih^ or ad>^ -h a^^ -|- Ozb^ 4- ad)^ -h ad>^. 

*=i 

Again we warn the reader that the superscript i in b^ is not the zth power 
of a variable b. 

The definition of the multiplication of two matrices can now be given 
in a neat form with the aid of the summation convention. Let 


al, 

’ * • , 

Ctm 

af, 

2 


a?, 

« 2 , ■ * ■ , 

ai 



K 

K 

7,2 ... 

7 

hi 

• • 

• « 

br 

■Lm 

7 02, • • • 

9mm 

m m m 

•L/m 

1 
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Then, by the product AB of the two matrices^ we shall mean the maZrix 


C - IM II , 

where 

4 = aihf (i « 1, 2, • • •, n; i = 1, 2, • • •, p). 

If cj is written out in extenso without the aid of the summation con- 
vention, we have 

4 = + • • • + ainhf. 


It should be emphasized here that, in order that the product of 
two matrices be well defined, the number of rows in the matrix B must 
be precisely equal to the number of columns in the matrix A . It follows 
in particular that, if A and B are square matrices of the same type, then 
AB as well as BA is always well defined. However, it must be empha- 
sized that in general J-B is not equal to BA, written as ^B 5^ BA, 
even if both AB and BA are well defined. In other words, matrix 
multiplication of matrices, unlike numerical multiplication, is not 
always commutative. 


Exercise 


The following example illustrates the non-commutativity of matrix 
multiplication. Take 


A = 


and 


Now 


B = 


Hence 


Similarly 


0 

1 


-1 

0 


1 

0 

0 

1 


so that a\ = 0, 0.2 — 1, a? = 1, = 0, 


so that h\ ~ —1, 62 = 0, = 0, = 1. 


cl = aii)? = (0)(-l) + (1)(0) = 0, 

4 = oifcS = (o)(o) + (i)(i) = 1 , 

cf = = (1)(-1) + (0)(0) = - 

4 - aStf = (1)(0) + (0)(1) = 0. 


AB - 


BA - 


1 , 


0 

1 

-1 

0 

0 - 

-1 

1 

0 


But obviously AB BA. 

The unit matrix I of order n has the interesting property that it 
commutes with all square matrices of the same order. In fact, if A is 
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an arbitrary square matrix of order n, then 

AI A. 


The multiplication of row and column matrices with the same 
number of elements is instructive. Let 


A = \\ ai II 

be the row matrix and 

B = II V II 

the column matrix. Then AB = aih\ a number, or a matrix with one 
element (the double-bar notation has been omitted). 


Exercise 


If ^ = II 1, 1, 0 II and B = 


0 

0 

1 


, then 


AB = (1)(0) + (1)(0) + (0)(1) = 0. 


This example also illustrates the fact that the 'product of t'wo 'matrices 
can he a zero matrix aUhoiigh neither of the rnultiplied matrices is a zero 
matrix. 

The multiplication of a square matrix uoith a column matrix occurs 
frequently in the applications. A system of n linear algebraic equations 
in n unknoums x^, x^, • • • , a:” 

ajx^ = h^ 

can he 'written as a single matrix equaiion 

AX = B 


in the unknown column matrix X = || || and the given square 

matrix A = || a} || and column matrix B \\ I]. 

A system of first-order differential equations 


dx^ 


= ajX^ 


can be written as one matric differential equation 


dX 

dt 


= AX. 


Finally a system of second-order differential equations occurring 
in the theory of small oscillations 

di^x^ V . 
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can be written as one matric second-order differential equation 



The above illustrations suffice to show the compactness and sim- 
plicity of matric equations when use is made of matrix multiplication. 


Exercises 

1. Compute the matrix AR when 


A = 

1,0 

0,1 

and B = 

0,0,2 


1,0 


-1,0,0 


Is Bi defined? Explain. 

2. Compute the matrix AX when 



1,3,0 


1 

A = 

-1, 2, 1 

and 1 = 

0 


0,0,2 


-2 


Is li defined? Explain. 
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ALGEBRAIC PRELIMINARIES (Continued) 

Inverse of a Idatrix and the Solution of Linear Equations.^ 

The inverse or reciprocal, of a real number a is well defined if 
a 0. There is an analogous operation for square matrices. If A is a 
square matrix 

^ = II II 

of order n and if the determinant | a} | 0, or in more extended no- 

tation 

^1 j * * * ? ^ 

a?, 

9^ 0 , 


(hr ' 'i (h \ 

then there exists a unique matrix, written A~^ in analogy to the inverse of 
a number, with the important properties 

(7 is the unit matrix.) 

The matrix A~^, if it exists, is called the inverse matrix of A, 

In fact, the following more extensive result holds good. A necessary 
and sufficierU condition that a matrix A - \\ a} || have an inverse is that 
the associated determinant a} \ 9 ^ 0 . 

From now on we shall refer to the determinant a = \ a} | as the 
determinant a of the matrix A, Occasionally we shall write | A | for 
the determinant of A, 

The general form of the inverse of a matrix can be given with the 
aid of a few results from the theory of determinants. Let a = | a) 
be a determinant, not necessarily different from zero. Let a} be the 
cofactor f of a{ in the determinant a; note that the indices i and j are 
interchanged in a} as compared with a{. Then the following results 

t The (n - l)-rowed determinant obtained from the determinant a by striking 
out the jih row and ith column in a, and then multiplying the result by (-1)*'^'^. 
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\a-^A =7 
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come from the properties of determinants : 

c^jOL% = a (expansion by elements of ^th row) ; 

o^-a% — a dk (expansion by elements of /bth column). 


If then the determinant a 9 ^ 0 , we obtain the following relations, 


( 2 * 2 ) 

on defining 


Vk = dl 




« • 


a* 


Let A = II a} ||, 5 = || 0} ||; then relations 2*2 state, in terms of 
matrix multiplication, that 

AB = J, BA ^ I. 


In other words, the matrix B is precisely the inverse matrix A~^ of A. 

To summarize, we have the following computational result: if the 
determinant a of a square matrix A = || a) || is different from zero, then 
the inverse matrix A~^ of A exists and is given by 

A-'- = II j8} II , 

where = — and a} is the cofactor of a\ in the determinant a of the 

€b 

matrix A. 

These results on the inverse of a matrix have a simple application to 
the solution of n non-homogeneous linear (algebraic) equations in n 
unknowns x^, x^, • • • , x^. Let the n equations be 

afjX^ = 

(the numbers a} are given and the n numbers are given). On de- 
fining the matrices 

A = \\ 4 II, X=\\ x‘ II, B = II 6‘ II, 


we can, as in the first chapter, write the n linear equations as one matric 
equation 


AX = B 


in the unknown column matrix X. If we now assume that the de- 
terminant a of the matrix A is not zero, the inverse matrix A~^ will 
exist and we shall have by matrix multiplication 

A-\AX) = A-^B. 

Since A~^A = 1 and IX = X, we obtain the solution 

X = A-^B 

of the equation AX = B. In other words, if is the cofactor of a\ in 
the determinant a of A, then x"^ = ocfb^/a is the solution of the system 
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of n eqvxiUons a)x^' = 6® under the condition a ^ 0. This is equivalent 
to Cramer’s rule^ for the solution of non-homogeneous linear equations 
as ratios of determinants. It is more explicit than Cramer’s rule in 
that the determinants in the numerator of the solution expressions are 
expanded in terms of the given right-hand sides 6^, • • • , 6"" of the 

linear equations. It is sometimes possible to solve the equations o^x^ = 
readily and obtain x^ = X}6^‘. The inverse matrix A~^ to A = || a) I| 
can then be read off by inspection — in fact, = || X} ||. 

Practical methods, including approximate methods, for the calcula- 
tion of the inverse (sometimes called reciprocal) of a matrix are given in 
Chapter IV of the book on matrices by Frazer, Duncan, and Collar. 
A method based on the Cayley-Hamilton theorem will be presented at 
the end of the chapter. 

A simple example on the inverse of a matrix would be instructive at 
this point. 

/ 

Exercise 


Consider the two-rowed matrix 

A = 


0 1 

-1 0 


According to our notations 

a\ = 0, aj = 1, a? = -1, 


Hence the cofactors otj of A will be 

«i = (cofactor of a\) =0, ocl — (cofactor of af) = — 1, 
Oil = (cofactor of al) = 1, cal = (cofactor of a^) = 0. 


Now A ^ = II 0j II , where /S} = a} /a. But the determinant of A is 
a = 1. This gives us immediately jSj = 0, = —1, jS? = 1, ^ = 0. 


In other words. 



Approximate numerical examples abound in the study of airplane- 
wing oscillations. For example, if 

0.0176, 0.000128, 0.00289 

A = 0.000128, 0.00000824, 0.0000413 , 

0.00289, 0.0000413, 0.000725 

then approximately 

170.9, 1,063., -741.7 

A-i = 1063., 176,500., -14,290. 

-741.7, -14,290., 5,150. 

3ee exercise 2 at the end of Chapter 7, 
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From the rule for the product of two determinants,® the following 
result is immediate on observing closely the definition of the product 
of two matrices; 

If A and B are two square matrices with determinants a and h respec- 
tively ^ then the determinant c of the matric product C - AB is given by the 
numerical multiplicaiion of the two numbers a and b, i.e., c = ah. 

This result enables us to calculate immediately the determinant of 
the inverse of a matrix. Since AA~^ = I, and since the determinant of 
the unit matrix I is 1, the above result shows that the determinant of 
A~^ is 1/a, where a is the determinant of A. 

From the associativity of the operation of multiplication of square 
matrices and the properties of inverses of matrices, the usual index 
laws for powers of numbers hold good for powers of matrices even 
though matric multiplication is not commutative. By the associativity 
of the operation of matric multiplication we mean that, ii A, B, C are 
any three square matrices of the same order, then t 

A(BC) = (AB)C. 

If then A is a square matrix, there is a unique matrix A A • • • A with 
s factors for any given positive integer s. We shall write this matrix 
as A® and call it the sth power of the matrix A. Now if we define 
A® = I, the unit matrix, then the following index laws hold for all 
positive integral and zero indices r and s: 

A^A* = A^A*- = A*+’' 

(AO* = (A*)’’ = A’’*. 

Furthermore, these index laws hold for all integral r and s, positive or 
negative, whenever A~^ exists. This is with the understanding that 
negative powers of matrices are defined, as positive powers of their inverses, 
i.e., A”'" is defined for any positive integer r by 

A-’- = (A-O'"- 

Multiplication of Matrices by Numbers, and Matric Polynomials. 

Besides the operations on matrices that have been discussed up to 
this section, there is still another one that is of great importance. If 
A = 11 a} 11 is a matrix, not necessarily a square matrix, and o: is a 
number, real or complex, then by aA we shall mean the matrix H cko} ||. 
This operation of multiplication by numbers enables us to consider 
matrix polynomials of type 

(2-3) aoA” + aiA”-^ + a2A«-2 + • • • + an-iA + aj. 

t Similarly, if the two square matrices A and B and the column matrix X have 
the same number of rows, then {AB)X = A{BX). 
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In expression 2 • 3, oro, ori, • • • , an are numbers, A is a square matrix, 
and I is the unit matrix of the same order as A. In a given matric 
polynomial, the aiS are given numbers, and A is a variable square 
matrix. 

Characteristic Equation of a Matrix and the Cayley-Hamilton Theorem. 

We are now in a position to discuss some results whose importance 
cannot be overestimated in the study of vibrations of all sorts (see 
Chapter 6). 

If A = II Gj II is a given square matrix of order n, one can form the 
matrix XZ — A, called the characteristic matrix of A. The determinant 
of this matrix, considered as a function of X, is a (numerical) poly- 
nomial of degree n in X, called the characteristic function of A. More 
explicitly, let /(X) = | XZ - A | ; then /(X) has the form /(X) =» X** + 
aiX"“^ + * * • + ctn-iX + Gn. Since an = /(O), we see that an = | —A | ; 
i.e.. On is (—1)” times the determinant of the matrix A. The algebraic 
equation of degree n for X. 

/(X) = 0 

is called th-e characteristic equation of the matrix A, and the roots of the 
equation are called the characteristic roots of A. 

We shall close this chapter with what is, perhaps, the most famous 
theorem in the algebra of matrices. 

The CayijEY-Hamilton Theorem. Let 

f(X) = X*' -j- aiX”~^ + • • * + On— iX + Gn 

he the characteristic function of a matrix A, and let I and O he the unit 
matrix and zero mnirix respectively with an order equal to that of A . 
Then the matric polynomial equation 

4“ ai.X^~^ H- * * • 4- Gn—xX 4- GnJ — O 
is saMsfied hy X — A. 


Example 


Take A = 


0 1 
1 0 


n = 2, and ai = 0, 02 = 


; then /(X) = 
-1. ButA2 = 


X -1 
-1 X 


1 

0 

0 

1 


= X2 - 1. Here 
Hence A^ —I =?; O. 


The Cayley-Hamilton theorem is often laconically stated in the 
form matrix satisfies its own characteristic equation.” In symbols, 
if /(X) is the characteristic function for a matrix A, then /(A) = O. 
Such statements are, of course, nonsensical if taken literally at their 
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CHARACTERISTIC EQUATION OF A MATRIX 

face value. However, such mnemonics are useful to those who thor- 
oughly understand the statement of the Cayley-Hamilton theorem. 

A knowledge of the characteristic function of a matrix enables one 
to compute the inverse of a matrix, if it exists, with the aid of the Cayley- 
Hamilton theorem. In fact, let A be an n-rowed square matrix with an 
inverse A~~^. This implies that the determinant a of A is not zero. 
Since Q ^ an = (— we find with the aid of the Cayley-Hamilton 

theorem that A satisfies the matric equation 

I 

I = - — -H + • • • -f- an-zA^ + an-iA]. 

an 

Multiplying both sides by A~^, we see that the inverse matrix A~^ can 
be computed by the following formula: 

(2*4) A-i = — -1- aiA^-2 + • - • -j- On_2A -l- an-xl']. 

dri 

To compute A~^ by formula 2*4 one has to know the coefficients ai. 
02 , • • a»-i, an in the characteristic function of the given matrix A. 

Let A = II Uj- II ; then the trace of the matrix A, written tr (A), is 
defined by tr (A) = a\, the sum of the n diagonal elements a}, • • •, 

o”. Define the numbers^ si, S 2 , • ' • , Sn by 

(2-5) si = tr(A), S 2 -tr(A 2 ), • • s,. = tr (AO, * * *, s« = tr (A”) 

so that Sr is the trace of the rth power of the given matrix A . It can be 
shown® by a long algebraic argument that the numbers ai, •••, o» 
can be computed successively by the following recurrence formulas: 

'ax = — Si 
02 = — + Sa) 

Os = —- 3(05281 05iS2 + S 3 ) 

(2-6) < ' 

1/ X 

an = — “(05n— iSi -H an—zS2 -!-••'+ 05iSn— 1 H" S^). 

L n 

We can summarize our results in the following rule for the calculation 
of the inverse matrix A~^ to a given matrix A. 

A Kule for Calculation of the Inverse Matrix A~^. First 
compute the first n — 1 powers A, A^, A””^ of the given n- 

rowed matrix A. Then compute the diagonal elements only of A**. 
Next compute the n numbers Si, S 2 , * * ’, Sn as defined in 2-5. Insert 
these values for the Si in formula 2*6, and calculate ax, o^, •**, an 
successively by means of 2-6. Finally by formula 2 • 4 one can calcu- 



14 


ALGEBRAIC PRELIMINARIES 


late A~^ from the knowledge of cti, • ■ a^, and the matrices A, A^, 

• • Notice that the whole is not needed in the calculation 

but merely Sn = tr {A^), the trace of A'^. 

Punched-card methods can be used to calculate the powers of the 
matrix A. The rest of the calculations are easily made by standard 
calculating machines. Hence one method of getting numerical solutions 
of a system of n linear equations in the n unknowns 

a}x^' = 6® (| a} I 0) 

is to compute of .4 = || a} || by the above rule with the aid of 
punched-card methods and then to compute A~^Bj where B = || 6* ||, 
by punched-card methods. The solution column matrix X = || || 

is given hy X = A~^B. 


Exercises 


1. Calculate the inverse matrix to -4 = 
Solution. 


0 1 
1 0 


by the last method of this chapter. 


A2 = 


1 

0 

0 

1 


, Si = 0, S2 = 2, ai « 0, ^2 -— 1 . 


Now A-^ =s — — \^A + ai/3 = A. Hence 

(Z2 


A-^ = 


0 1 
1 0 


2 . See the exercise given in M. D. Bingham's paper. See the bibliography. 

3. Calculate A~^ by the above rule when 


15 

11 

6 

-9 

-15 

1 

3 

9 

-3 

-8 

7 

6 

6 

-3 

-11 

7 

7 

5 

-3 

-11 

17 

12 

5 

-10 

-16 


After calculating A^, A®, A*, and the diagonal elements of A®, calculate Si = 6 , 
S 2 = —41, S 3 = —217, S 4 = —17, Ss = 3185. Inserting these values in 2-6, find 

Oi “ — 5, Qi2 = 33, tt 3 = — filj ^4 = 135, ds — 225. 


Incidentally the characteristic equation of A is 

/(X) = X® - 5X* + 33X3 _ 51x2 + 135X + 225 
= (X + 1)(X2 - 3X + 15)2 * 0. 

Finally, using formula 2-4, find 

-207 64 -124 111 171 

-315 30 195 -180 270 

-315 30 -30 45 270 

-225 75 -75 0 225 

-414 53 52 -3 342 




CHAPTER 3 

DIFFERENTIAL AND INTEGRAL CALCULUS OF 

MATRICES 


Power Series in Matrices. 

Before we discuss the subject of power series, it is convenient to 
make a few introductory remarks on general series in matrices. Let 
Ao, Ai, A2, As ■ * ' be an infinite sequence of matrices of the same type 
(i.e., same number of rows and columns) and let = Ao + Ai + A2 
-I- • • • + Ap be the matric sum of the matrices Aq, Ai, A2, • • • , and Ap. 
If every element in the matrix Sp converges (in the ordinary numerical 
sense) as p tends to infinity, then by /S = lim ^Sp we shall mean the 

p— >00 

matrix S of the limiting elements. If then the matrix S = lim Sp exists 

p — >00 

in the above sense, we shall say, by definition, that the matric infinite series 

00 

^A, converges to the matrix S. 

r=0 

Example 

1 1 ^ 1 

Take Ao = I, Ai = J, A2 = ~J, A3 =* ~I, • • - , Ai =» —J, * * * . Then 

juI o I ^ I 

iSp = Ao 4- Ai + A 2 + • • • + Ap “ "b 1 "b + ■ * * + 

Hence, on recalling the expansion for the exponential e, we find that 

CO 

lim Sp = el. In other words, Ea. = el. 

p — > CD r— 0 

If A is a square matrix and the ai, 02, • • * are numbers, one can 

consider matric power series in A 

00 

r=0 

In other words, matric power series are particular matric series in which 
each matrix At is of special type \ Ar = arA^, where A*' is the rth power 
of a square matrix A. (A® = / is the identity matrix.) Clearly matric 
polynomials (see Chapter 2) are special matric power series in which all 
the numbers at after a certain value of i are zero. 

An important example of a matric power series is the matric expo- 
nential Junction e^ defined by the following matric power series: 

4 +••• + •■ •• 


t The index r is not summed. 


15 
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The following properties of the matrix exponential have been used 
frequently in investigations on the matrix calculus: 

1. The matric power series expansion for is convergent^ for all 
square matrices A. 

2. whenever A and JB are commutative matrices, 
i.e., whenever AB = BA. 

3. e^e ^ — e~^e^ = 1. (These relations express the fact that 
is the inverse matrix of e^.) 

Every numerical power series has its matric analogue. However, 
the corresponding matric power series have more complicated proper- 
ties for example, . Other examples are, say, the matric sine, 

A, and the matric cosine, cos Aj defined by 

sin A = A — —A® + — . 

3! 5! 

cos A = 7 - ^A^ + \-A^ - • • • . 

2! 4! 

The usual trigonometric identities are not always satisfied by sin A 
and cos A for arbitrary matrices. 


Differentiation and Integration of Matrices Depending on a Numeri- 
cal Variable. 

Let A{t) be a matrix depending on a numerical variable t so that 
the elements of A (t) are numerical functions of t. 


A(f) 


a\{t), ai{t) 

alit) 


II <(0, <(0 

Then we define the derivative of A(t)j and write it hy 

dt 


da\(t) dc^jt) da^{t) 

dt ^ dt ’ ^ dt 

da\(t) dd^{t) dai{t) 

dA (0 dt ' dt ' ' dt 

dt 


dafCO da^(t) da^(t) 

dt ' dt dt 
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Similarly we define the integral of A(t) by 

SaXit) dt, dt, 

Salit) dt, Sc^it) dt, 


S A (t) dt 


, SoL^Sfi) dt 

•y dt 


y*ar(0 dt, dt, • • S an {t) dt 

It is no mathematical feat to show that differentiation of matrices 
has the following properties: 


(3-1) 

(3-2) 

(3-3) 

etc. 


d\iA (0 + B(0I1 dA (t) ^ dB(t) 


dt dt 

dlAit)B(t)'} dA(t) 


dt 


dt 


dt 


Bit) + A(t) 


dB(f) 

dt 




+ A(e)jS(0 


dCjt) 
dt ^ 


There are important immediate consequences of properties 3 • 2 and 
3-3. For example, from 3-2 and A~^{t)A(t) = I, we see that 

dA-^(t) 


(3-4) 


dt 


-A-^(t)^^^A-Ht). 


dt 


it). 


Also, from 3 • 3, we obtain 
dA^it) dA(t) 


(3*5) 


dt 


dt 


AKt) + A(fr:^A(t) + 


dt 


dt 


There are similar formulas for the derivative of any positive integral 
power of A (t) . 

If i is a real variable and A a constant square matrix, then one obtains 

d(rA) 


dt 


rf’^^A. 


Then, with the usual term-by-term differentiation of the numerical 
exponential, the following differentiation can be justified: 


(3-6) 


d (c^^) t^ t^ 


+ • • * 


== Ae*'^ = d^A. 


tA 
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There is an important theorem in the matrix calculus that turns up 
in the mathematical theory of aircraft jivUer (see Chapter 7). The 
proof, into which we can not enter here, makes use of the modem 
theory of functionals. 


Theorem. If FQC) is a power series that converges for all X, then the 
maJtric power series F(^A) can he computed hy the expansion^ 

(3-7) F(A) = Y,F(}.i)Gi, 


where A is an n-rowed square matrix with n distinct characteristic roots 
Ai, Az, • • ‘ , A„, and Gx, 0%, • • ’ j Gn are n matrices defined hy^ 


(3-8) 



jjiii 


— 


There are a few matters that must be kept in mind in order to have 
a clear understanding of the meaning of this result. In the first place 
the matric power series F(A) == <xoI + axA + ol^A^ + • * * + * • ■ when- 
ever F(A) = ao + <x{K -h -!-••• + ••'• In other words A® = 1 is 
^‘replaced” by A^ = 7, the unit matrix, in the transition from 7"(A) to 
F(jA). Secondly to avoid ambiguities we must write explicitly the 
compact products occurring in equation 3 • 8. 


I I (At — Ai) = (Ai — Ai)(A2 — Af) * ■ * (Ai— 1 — Ai)(Af+i A*) • • • (An Ai), 


5 


(Ay7 - A) = (Ai7 - A)(A27 - A) • - - (A»-i7 - A){\i+xl - A) • • ■ 

(An7- 


A). 


There are special cases of particular interest in vibration theory (see 
Chapters 6 and 7). They correspond to the power of a matrix A^ and 
the matrix exponential e^. The expansion 3-7 yields immediately 


(3-9) A’- = £x;(?j 

»=«1 

and 

(3-10) = 'f^e^Gi, 

i=\ 

where the matrices G» have the same meaning as in 3 • 8. 


Exercise 


Calculate the matrix when A is the matrix A 
by calculating directly. 


0 

1 

1 

0 


Chock tho result 
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„ The ckacteristic roots aie \i = 1, Aj = -1. Hence the matrices 


III and (?2 arc as follows: 


W-i 1 


1 




U-i 1 


1 




1 1 
1 1 
1 -1 
-1 1 






e 


a » 


2 


1 

+- 

1 -1 

h 1 


-1 1 


I 



coshl sinhl 
slnhl coshl 



CHAPTER 4 


DIFFERENTIAL AND INTEGRAL CALCULUS OF 
MATRICES (Continued) 


Systems of Linear Differential Equations with Constant Coefficients. 


The matric exponential has important applications to the solution of 
systems of n linear differential equations in n unknown functions 

and with constant coefficients a}. The variable t 
is usually the time in physical and engineering problems. Without 


defining the derivative 


dt' 


we merely mentioned in the first chapter 


that we can write such a system of equations as one matric equation 


(4.1) 


dX{t) 

dt 


= AX{t). 


Having defined the matric derivative, we are enabled to view this 
equation with complete understanding. 

From formula 3 • 6 of the previous chapter we find that 

( 4 . 2 ) 

dt 


where U is an arbitrarily given value of t. But this result is equivalent 
to saying that X{t) = is a solution of the matric differential 

equation 4 • 1 for an arbitrary column matrix Xq. A glance at the 
expansion for the matric exponential shows that the solution 

X(0 has the property 

X(h) = Xo. 

In summary, we have the result ^ that 
(4*3) X(0 = 

is a solution^ o/ 4 • 1 with the 'property that Xit^) = Xq for any preassigned 
constant column moirix Xq. 


so that 


Example 


dx^it) 

dt 



dx\t) 

dt 


= x^ 


dt 


AX, 
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where 


A = 


0 1 
1 0 


and X = 




Now Xi = 1, Xg = — 1, 
chapter that 



and we saw in the last exercise of the previous 


1 1 
1 1 



-1 

1 


Hence 


Z=1 


cosh {t — io) sinh {t — ’ 

sinh {t — <o) cosh (t — io) 


Therefore the unique solution of the differential system 

dt 


= AX, X(to) = Xo 


xl 

xl 


is 


Xit) = 


cosh {t — ^o) sinh {t — to') 
sinh (f — io) cosh (t — to) 


This means that the unique solution of the differential system 
dx^ . dx^ 


dt 


x^. 


dt 


= xK 


x^(to) = xi, x^(to) = x% 


= [[cosh — ^o)!]3^o + Csinh {t — <o)Ila:§ 
x^lt) = [[sinh (t — io)3iCo + [[cosh {t — to)~]x%. 


Systems of Linear Differential Equations with Variable Coefficients. 

Although the matric exponential is not applicable to the solution of 
a system of linear differential equations with variable coefficients 
a)it), there are some analogous matric expansions that enter into the 
solution of such a system. The system of differential equations 

dxKt) 

(4-4) = a){t)xKt) 


is written as one matric differential equation 


(4*5) 


dX(t) 

dt 


= A(0X«) 


where A{t) — || |1 and X(f) is the column matrix of the n un- 

known functions x'^(t). 

On integrating both sides of 4*5 between U and t we obtain the 
equivalent matric equation 

(4-6) X(t) = X(to) + 


r A (s)X(s) ds. 

Jta 
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By the method of successive substitutions, we are led to consider the 
following expansion as a solution of 4*6: 


(4-7) X(t) = [J + r^(s) ds + CaCs) ds T^Csi) dsx 

Uio *Jto 


+ 


+ • • •] 

X(^o). 


NTow the method of successive substitutions for equation 4-6 can be 
described as follows. In the integral term in 4 • 6 substitute for X(s) its 
equivalent as given by formula 4 • 6 itself. This yields 


X(ty^ X(to) + r CaCs) ds~\x(to) + r^(s) ds rA(si)XCsi) dsi. 

— .V fo J to to 


Again substituting for X(si) its equal as given by 4-6 we are led to a 
new expansion for X(0. Continuing indefinitely this way we are led 
to the matric infinite series 4 • 7. 

If we define the matrix 


(4-8) iliiA) = 7 4- r A(s) ds + f A(s) ds TaCsi) dsi 

•J to d to d to 

J ft pa psi 

A (s) ds I A (si) dsi I A (§ 2 ) dsz + * * * + • • • , 

then it can be proved that, for a}(0 continuous in i^o ^ ^ ^ ^ 1 , 

(4-9) X(0 = i4(A)Xo 

is the unique solution of the matric differential equaiion 4-5 that takes 
on the arbitrarily given constant matric value Xq for t *= ^o- It is often 
simpler to carry out the matrix multiplications first in 4 • 8 and 4 • 9 
before carrying out the successive integrations. If the matrix is inde- 
pendent of t, then, by an evident calculation, solution 4-9 reduces 
precisely to the matrix exponential type 4-3. 

For approximate numerical calculations, a few terms in the expansion 
for Qfti(A) may suflBce in 4-9 to give a good approximation to the 
solution of the matric differential equation 4-5. 


Exercises 


1. Integrate by matrix methods the second-order differential equation 

cPxit) 




— x(t) >= 0 


subject to the initial conditions x(jto') - Xo, = yo. 

^97rite the differential equation as a system of two first-order equations 


dx^ 

dt 


x^ 
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with initial conditions x^{U) « Xq, x\t^ - 2 / 0 , and use the results of the example il- 
lustrating formula 4*3.) 


2. Integrate by matrix methods the second-order differential system for har- 
monic oscillations with frequency W 


dhjt) 

If 


+ = 0 , 




= 2 / 0 . 


(JjTmi. Write the equation as a system of two first-order linear equations.) 

3 . Discuss the solutions of the differential equation 


dH dx , 


0 


for free damped oscillations by matrix methods with the restriction that jS 9 ^ 2\/ km. 
m = mass, jS = damping factor, and k =* elastic constant, so that all three j3, k 
are positive constants. Clearly the restriction rules out the critical damping case. 
{EinL Write the differential equation as a first-order matric differential equation 



where 



0 

1 


k 



m 

m 


and notice that the characteristic equation of this matrix is the “characteristic 
equation” of the given second-order differential equation in the usual elementary 
sense.) 

4 . Integrate by matrix methods the second-order differential equation 


d^x{i) 

df 


— tx{t) = 0 


subject to the initial conditions a;(0) « xqj 



= 2/0. 



CHAPTER 5 


MATRIX METHODS IN PROBLEMS OF SMALL 

OSCILLATIONS 


Differential Equations of Motion.^ 

The problem of small oscillations^ (of conservative dynamical sys- 
tems) about an equilibrium position concerns itself with the solution 
of the Lagrangian differential equations of motion in which the kinetic, 
and potential energies are homogeneous quadratic forms, in the veloci- 
ties and coordinates respectively, with constant coefficients. The 
theory is approximate in that the constancy of the coefficients in the 
kinetic energy and the quadratic type of the potential energy are due 
to approximations in the actual form of the kinetic and potential 
energies respectively. If, without loss of generality, we take all the 
coordinates of the equilibrium position to be zero, these approximations 
are due to the assumed smallness of the coordinates and velocities 
about the . equilibrium position. 

Let 


and 


_ - dq* 

y = Q^ij = M 


be the kinetic and potential energies respectively of our oscillating 
system with n degrees of freedom. In view of what we have already 
said, the Uy and 6^ are constants. We shall consider the case in which 
the equilibrium point is stable, i.e., the potential energy V has a mini- 
mum at g* = 0. Now it can be proved that the 'positive definiteness of V 
is a necessary and sufficient condition that (0, 0, • • • , 0) be a stable 
equilibrium point. V is, by definition, positive definite if F > 0 for 
all g' and F = 0 if and only if = 0. Clearly the kinetic energy T 

dq^ 

is positive definite in the velocities — 

dt 

Lagrange’s equations of motion for our oscillating system are 


dt\bqy 




a == 1, 2, • . n) 


on using the notation 



If we use the explicit form for the 
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kinetic and potential energies, Lagrange^s equations reduce to the 
system of n second-order differential equations 

(5*1) 


where we have used the notation q^' = 
matrices f 

A == II a^j II , 
and the unknown colunm matrix 


dt^ 


If we define two square 


B = II bij 


Qit) = II qKt) II , 

then we can write our differential equations 5 * 1 of motion as the one 
matric differential equation 


(5-2) 


dt^ 




Since the kinetic energy is positive definite, it can be proved that 
the determinant | A \ 9 ^ 0. Hence it follows from our discussion in 
Chapter 2 that the inverse matrix A~^ exists. On multiplying both 
sides of equation 5 • 2 on the left by A~^ and remembering that A~^A — I, 
the unit matrix, we obtain the following equivalent matric differential 
equation 


(5-3) 


d'^Qjt) 

df- 


-CQ(t), 


>■ 


where C is the (constant) square matrix C = A~^B. 

To summarize, we have the following result. If A and B are the 
constant square matrices of the coefficients of the kinetic and 'potential 
energies respectively ^ then the motion of our oscillatory system is governed 
by the rrwiric differential equation 5-3. 


Illustrative Example 

Two equal masses, each of mass m, are connected by a spring with 
elastic constant h while each mass is connected to a fixed wall by a 
spring with elastic constant k. The kinetic and potential energies of 
this two-degree-of-freedom problem are 

^ - f Kf )■ - (f )■] 

V = + (g* - g^)^. 


t Recall the notations for matrices given in Chapter 1. 
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where and are the respective displacements of the centers of the 
two masses ‘‘ parallel” to the springs and are measured from the 
equilibrium position in which all three springs are unstressed. Hence 
by a direct calculation from the kinetic and potential energies we find 
that the matric equations of motion are 

(5-4) = -CQ(t), 

where 

il 2k k II 


= -CQ(«), 


C = 


dO(t') 

If we define the column matrix i2(0 == 2 / ^ "”3 — ~ we 

T \Jij Cut 

can write the second-order matric differential equation 5*4 as a first- 
order matric differential equation 

(5-5) = US(t), 


(5-5) 


= US{t), 


where 


^(0 = 


qKt) 

QKt) 

r\t) 

r\t) 


and 


(5-6) 


U = 


0 

0 

-2k 


— 00 


k -2k 


m m 


0 0 


The characteristic equation of the matrix XJ turns out to be 

X4 + — + — = 0. 


Now — > 0, so that there are four distinct pure imaginary character- 
istic roots of U given by 


(5-7) 


-Zk 


, X 2 = — 


-Zk 


i X3 = 


, X4 = — 



27 


Eierdse 


A slialt of leDgth 2f, kd at one end, carries one disk at the free end and another^ 
in the middle. If ^ is the moment of inertia of each disk, and f, f are the respec- 
tive angular defections of the two disks, then the kinetic and potential energies are 





under the assumption that the shaft has a uniform torsional stiffness T. Find the 
matric differential equation of motion. Write this equation as a system of two 
firstK)rder matric equations. Discuss the solutions of this system and then the mo- 
tion of the disks. 



CHAPTER 6 


MATRIX METHODS IN PROBLEMS OF SMALL 
OSCHLATIONS (Continued) 


Calculation of Frequencies and Amplitudes. 

Let us inquire into the pure harmonic solutions of 
equation of motion 


( 6 * 1 ) 


d^Q(V) 

dt^ 


-CQ(t\ 


our differential 


where C = We thus seek solutions of 6-1 of type 

(6*2) Q(t) = sin (cot + ^)r, 

where co is an angular frequency, an arbitrary phase angle, and V 
a column matrix of amplitudes. On substituting 6*2 in 6* 1 we obtain 

— co^ sin (cot + ^)r = —sin (cot + ^)Cr. 

Hence a necessary and sufficient condition that 6*2 be a solution of 
6 ■ 1 is that the frequency co and the correspon ding column matrix T of 
amplitudes satisfy the matrix equation 

(6*3) (coU-C)T=^0. 

In order that there exist a solution matrix P 5 *^ 0 of 6*3, it is clear 
from the theory of systems of linear homogeneous algebraic equations 
that 0)2 must be a characteristic root of the matrix C. Since C = A-^B, 
we verify immediately the statement that 

(6-4) (oU -C = A-^(co^A - B). 

On recalling that the determinant of the product of two matrices is 
equal to the product of their determinants, we see that the determinant 

I - B) I = I 4-1 I 1 - B I 

and hence, by 6*4, 

I 0)2/ - c I = I A-1 I I 0)2A - R I . 

But 1^ A~^ I ?£ 0, so that the characteristic roots of the matrix C are 
identical with the roots of the frequency equation 

(6-5) I XA-R 1 =0. 

Since the kinetic and potential energies are positive definite quadratic 
forms, it can he proved (see any book on dynamics such as Whittaker's) 
that all the roots of the frequency equation are positive. Hence all the 
characteristic roots of the matrix C — A~^B are positive. 
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Since the potential energy is a positive definite quadratic form, it 
follows that I J5 I > 0 and hence that exists. Therefore C~^ 
exists and is given by C~^ = B~^A. On multiplying both sides of equa- 
tion 6*1 by C~^, we obtain the matric differential equation 


(6-6) 


Q{t) 


-Z) 


d^Q{t) 

dt^ 


where D = = B-^A. This equation is obviously equivalent to equa- 

tion 6T. If we now proceed with 6*6 as we did with 6*1, we are led 
to the equation 

(6-7) (i/-i>)r-o. 

This equation can also be derived by operating directly with equation 
6-3. It is clear from equations 6-3 and 6-7 that the characteristic 
roots of the matrix D = are the corresponding reciprocals of the 
characteristic roots of the matrix C. We shall call I) the dynamiad 
matrix. 

Let us write 6*7 in the equivalent form 
(6-8) r = co^DT. 

The classical method of finding the frequencies and amplitudes of 
our oscillating system consists in first finding the frequencies by solving 
the frequency equation 6 ■ 5, or equivalently in finding the characteristic 
roots of the matrix C = A~^B, and then in determining the ampli- 
tudes by solving the system of linear homogeneous equations that 
corresponds to the matrix equation 6-3. Such a direct way of 
calculating the frequencies and amplitudes often involves laborious 
calculations. For approximate numerical calculations, the method 
of successive approximations when applied to equation 6-8 greatly re- 
duces the laborious calculations. This is especially true when only the 
fundamental frequency (lowest frequency) and the corresponding 
amplitudes are desired. ^ We shall assume now that all the frequencies of 
our oscillating system are distinct. The method of successive approxi- 
mations for equation 6-8 is as follows. Let To be an arbitrarily given 
column matrix. Define 

Fi = ay^DVo 

Fa == oj^DFi 


and in general 


Tr = co'^DVr-i. 
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By successive use of this recurrence relation we can express Fr in terms 
of To. In fact, we have 

(6*9) r,. = (co2)»-D’To, 

where is the rth power of the dynamical matrix D. Now it can be 
shown that for large t the vatio of the eleTneyits of the coIultyiti vnottvix 
0 to the correspoTiding elements of the column matrix D^Tq is approxi- 
mately a constant eqztal to cof, the square of the fundamental frequency coi 
of our fundaTYiental mode of oscillation wzth distinct frequencies. The 
matrix To is only restricted by the nour-vanishing of R (see equality 6-11 
below). 

The proof of this result is a little involved and makes use of the 
Cayley-Hamilton theorem, a theorem "j* of Sylvester, and a few other 
theorems on matrices.^ These theorems are instrumental in showing 
that, for r large enough, the following approximate equality holds: 


( 6 - 10 ) 






(X2 — Xi)(X3 — Xi) • • • (Xn ~ Xi) 




oc 


where (a) Xi > X 2 > • • • > X» are the characteristic roots of Z>, and 
hence Xi = -^ in terms of the fundamental frequency coi; (fo) the num- 

COi 

hers O'!, • • - , are proportional to the amplitudes of the funda- 

mental mode of oscillation; and (c) 


( 6 . 11 ) 



In 6" 11, the To are the elements of the arbitrarily chosen column matrix 
To and the ^ 1 , • • *, are n constants that are themselves obtainable 
by a successive approximation method. 

Clearly, DTo is a column matrix. Hence the approximate formula 
6*10 shows that for r large enough, and for arbitrarily chosen To, such 
that R 5-5 0, the column matrix Z)»To has elements proportional to the 
amplitudes of the fundamental mode of oscillation. All this is subject to 
the restriction that all the frequencies are distinct. 


On using equation 6-3 instead of 6-7, one can similarly obtain the 
n ri(Xs/ — A) 

t P(A) - £^P(Xr)6fr, where Gr = and Xi, • . are characteristic 


s^r 


roots of A. See the discussion in Chapter 3. 
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greatest frequency and corresponding amplitudes of our oscillating 
system. The intermediate overtones and corresponding amplitudes 
can be obtained by the above successive approximation methods 
on reducing the number of degrees of freedom successively by one. 
Any one interested in these topics wiE find the following paper by 
W. J. Duncan and A. R. CoEar very useful: “A Method for the 
Solution of Oscillation Problems by Matrices,” Pfihopfiicd Ma^m 



CHAPTER 7 


MATRIX METHODS IN THE MATHEMATICAL THEORY 

OF AIRCRAFT FLUTTER 


In recent years a group of phenomena known under the caption 
“flutter has engaged the attention of aeronautical engineers. The 
vibrations taking place in flutter phenomena can often lead to loss of 
control or even to structural failure in such aircraft parts as wing, 
aileron, and tail. Such dangerous situations may arise when the 
airplane is flown at a high speed. It is of the greatest practical im- 
portance therefore so to design the plane as to have the maximum 
operating speed less than the critical speed at which flutter occurs. 
Unfortunately experiments in wind tunnels are idealized and difficult, 
and actual flight testing is obviously highly dangerous. It is here 
that mathematics enters the stage at a most opportune moment. 
Although the results of mathematical theories of flutter are now being 
applied in the design of aircraft, the need for an adequate mathe- 
matical theory is becoming critical. There is no time in this brief set 
of mathematical lectures to deal adequately with the present simplified 
mathematical theories of the mechanism of flutter. We shall only 
give the matric form for the equations of motion and say a few words 
about the approximate solutions with the aid of matrix iteration 
methods. 

The vibrations of an airplane wing and aileron can be considered as 
those of a mechanical system with three degrees of freedom : the bend- 
ing and twisting of the wing accounts for two degrees of freedom, and 
the relative deflection of the aileron gives rise to the third degree of 
freedom. Not only is the system non-conservative, but there is the 
additional complication of damping forces leading to terms depending 
on the velocities in the equations of motion. The differential equations 
of motion are of type 


(7*1) 




-h = 0, 


a system of three linear differential equations in the three unknowns 

Since there are three degrees of freedom, all indices 
have the range 1 to 3. The constant coefficients hij, Cij are com- 
puted from a large number of aerodynamic constants of our aircraft 
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structure; see T. Theodorsen^s ‘‘General Theory of Aerodynamic 
Instability and the Mechanism of Flutter,” N. A. C. A. Report 496, 
for 1934, pp. 413-433, especially pp. 419-420. Now the general 
structure of equations 7 • 1 differs from the equations of motion of the 
preceding two chapters in that 6^- 5^ bji (giving rise to a non-conserva- 

tive system) mid in the presence of the linear damping terms 

If we define A = \\ aij |I, jB = || 6^- ||, C = || Cij ||, Q(t) = || qKt) ||, 
then the equations of motion 7*1 can be written as the one matric 
differential equation 


(7-2) 


d^Q(i) dQ(f) 


dt^ 


dt 


in terms of the three known constant matrices A, S, (7, and the un- 
known column matrix Q(t). As we are interested in small oscillations 
around an unstable point of equilibrium, it is to be expected that 
complex imaginary frequencies will play a role in the work. 

Since A arises from the kinetic energy, A~^ exists and hence 7.2 is 
equivalent to 

d^Q if) ^ dQ({) 

de ~ ^ dt 


(7-3) 


- A-^BQ(t). 


We can replace the one second-order differential equation 7-3 by an 
equivalent pair of two first-order differential equations with the column 
matrices Q(t) and R(t) as unknowns 

) dt 
\ dt 




A-^BQ(t) - A-^CR(t). 


Define the column matrix of six elements 

Q(t) 


S(t) = 


Rit) 


and the constant square matrix of six rows 

O, / 

(7-6) 

I -A-^B, ~A-^C 

where O and I are the three-rowed zero and unit matrices respectively. 
Then equations 7 • 4 can be written as the one first.-order matric differen- 
tial equation 

dS(t) 


(7-6) 


dt 


= US(J0. 
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We are led therefore to consider solutions 

JS(t) = e^A (A, a constant column matrix of six elements) 
of 7 • 6. This obviously leads us to the equation 
(7-7) (X/-C7)A = 0, 

where I is the six-row unit matrix. To solve the problem we must 
get good approximations to the values of \ and the matrix A that will 
satisfy 7-7. The matrix iteration method for small oscillations of 
conservative systems can now be applied with some modifications made 
necessary by the fact that the possible values of X in 7 • 7 are in general 
complex imaginary. If Xi, X 2 , • • X„ are the characteristic roots of 
the matrix U lexicographically arranged so that their moduli f are in 
descending order, i.e., | Xi | > | X 2 | > • • • > | Xn | , then the 
characteristic root Xi with the largest modulus can be obtained by the 
methods of the previous chapter. Some further aids in computation 
of the real and imaginary parts of complex characteristic roots are 
given on pp. 148-150 and 327-331 of the Frazer-Duncan-Collar book 
on matrices. A more readable and self-contained account is given in 
the paper by W. J. Duncan and A. R. Collar entitled “Matrices 
Applied to the Motions of Damped Systems,’^ Phil. Mag., vol. 19 
(1935), pp. 197—219. An illuminating discussion of a specialized 
flutter problem with two degrees of freedom is given in a 1940 book by 
Kdrmdn and Biot, Mathematical Methods in Engineering, pp. 220—228. 
It would be interesting and instructive to solve such specialized flutter 
problems with the aid of the matrix calculus. 

Another useful method of solving flutter problems is the combi- 
nation of matrix methods and Laplace transform methods. The 
Laplace transform of a function x(f) is a function x(j>) defined by 

J r»co 

e~^*x{t) dl. 
o 

If one is willing to omit the proofs of one or two theorems, the whole 
Laplace transform theory needed does not require one to be conversant 
with the residue theory of complex variable theory. For such an 
elementary treatment of Laplace transforms see Operational Calculus in 
Applied Mathematics by Carslaw and Jaeger, Chapters I— III. The 
methods given there can be immediately extended in the obvious way 
to apply directly to the matric differential equations 7 • 3 for flutter 
problems. A good table of Laplace transforms together with mechan- 
ical or electric methods can cut down the labor of flutter calculations 

t The modulus o f a comp lex number z = x y is denoted by \z\ and is 

defined by j 2 ; | = + y^. 
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materially. Unfortunately there is no time to take up these matters 
in detail in our brief introductory treatment. 


Exercises 


1 . Show that the matric differential equation 7*6 for flutter can be written as 

S(t) == U- i dS(t) , 
dt 

where U~\ the inverse matrix of U, is given by the six-row square matrix 

-B-iC -B-^A 
I O 


u-^= 


2 . A model airplane wing is placed at a small angle of incidence in a uniform air 
stream. The three degrees of freedom are the wing bending, wing twist, and aileron 
angle measured relative to the wing chord at the wing tip. When the wind speed 
is 12 feet per second, the matrices for the differential equations of flutter are as 
follows. The data are obtained from R. A. Frazer and W. J. Duncan, “The Flutter 
of Aeroplane Wings,” Reports and Memoranda of the Aeronautical Research Com- 
mittee, No. 1165, August, 1928. 



17.6 

0.128 

2.89 

A = 

0.128 

0.00824125 0.0413 


2.89 

0.0413 

0.725 11 


121.042 

1.89 

15.9497 

B = 

0 

0.027 

0.0145 


11.9097 

0.364 

15.4722 


7.65833 

0.245 

2.10 

C = 

0.023 

0.0104 

0.0223 


0.60 

0.0756 

0.658333 


= matrix of damping coefficients. 


Show that 



0.170883 

1.06301 

-0.741731 

A-i = 

1.06301 

176.433 

-14.2880 


-0.741731 

-14.2880 

5.14994 1 


and that the matric differential equation for flutter is 


dSC) 

dt 


= US{t), 


where the “flutter matrix” U is given by 


0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

-11.8502 

-0.08168 

8.73520 

-0.888089 

0.003153 

0.105747 

41.4969 

-1.57195 

201.554 

-3.62604 

-1.01517 

3.23949 

28.4464 

-0.086931 

,-67.6433 

2.91908 

-0.059016 

-1.51412 1 


For the lengthy details of the calculations of flutter frequencies and amplitudes 
see the Phil. Mag. 1935 paper by Duncan and Collar. 
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More recent developments in aircraft design require an extension of 
the flutter theory to handle four-degree- (or more) of-freedom problems 
in which the motion of the tah defines the fourth degree t of freedom 
with generalized coordinate The aerodynamic forces and moments 
are obtained theoretically in accordance with T. Theodorsen’s and 
I. E. Garrick’s investigations and not from wind-tunnel data. It is 
for this reason that the coefiBlcients in the differential equations of 
motion will be in general complex. See Fig. 7*1. The four differential 
equations of motion are of the form % 


(7-8) 





where all the indices, in contradistinction to 7-1, have the range 1 to 
4 and the coefficients hij, and Cij are in general complex. 

The flutter velocity v appears in general in the coefficients c*y and 6yy 


and is replaced by the quantity — , where h, co, and k are the airfoil 

semi-chord, flutter frequency, and the flutter parameter A; = — respec- 

tively. This yields a system of four linear differential equations 7*8 

t In accordance with the flutter notation used in this country, k, = <x, 
g® = |8, g* — -y- See Fig. 7 • 1. 

t The differential equations of motion and the contributions of the aerodynamic 
forces and moments to the coefficients of the differential equations for the four- 
degree-of-freedom problem are given in the Douglas reports. 
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in. which Cij and are expressed as functions of the flutter frequency ca 
and not of the flutter velocity v. If one then considers solutions of the 
form 

2 ’' = {i = 

of the differential equations, one is led to a system of four linear alge- 
braic equations with complex coefl&cients, some of which are functions 
of CO and the structural damping coefl&cients gTy. The damping coeffi- 
cients Qij are defined by 

Fi = V-l QijKijq^' (gij = 0 if i^j), 

where Ft is the damping force in the g*th direction and the Kij are the 
spring constants. Upon obtaining the characteristic roots of the 
matrix of the linear algebraic equations by matric iteration methods, 
one ultimately finds the flutter velocity v as a function of the flutter 
frequency <a and the structural damping coefficients Qij. Flutter is 
likely to occur if the structural damping coefficients Qij from the 
pure imaginary part of the characteristic root exceeds 0.03, provided 
that no extraneous damping devices are used. The algebraic equa- 
tions are so arranged — and this constitutes an important aspect of 
the development in that it lends itself to matric iteration procedures — 
that the characteristic roots are of the form 

^--, + ig (* = vCI), 

where c is some constant, co is the flutter frequency, and gr is a structural 
damping coefficient. 

The flutter analyst is attempting to approximate the actual flutter 
characteristics of the airplane by representing them by as small a 
number of degrees of freedom as possible. The design of faster and 
larger aircraft requires the consideration of a larger number of degrees 
of freedom so as to make the flutter analysis an adeqvMe approximation. 
When many degrees of freedom are required to represent the flutter 
characteristics of an airplane, the need for matrix methods becomes 
acute. Matrix methods t also serve to improve the theory of the 
mechanism of flutter. 

t Matrix methods are also used in treating other phenomena related to flutter. 
See Douglas reports. 



CHAPTER 8 

MATRIX METHODS IN ELASTIC DEFORMATION 

THEORY 

Although the tensor calculus is the most natural and powerful mathe- 
matical method in the treatment of the fundamentals of elastic def- 
ormation of bodies, the matrix calculus can also be used to advantage 
in furnishing a short and neat treatment. This chapter is purely 
introductory and suggestive. 

Let a medium be acted on by deforming forces. The position of the 
medium before and after deformation will be called the initial and 
final state of the medium respectively. Let a? be the rectangular 
cartesian coordinates of a representative particle of the medium in the 
initial state, and the rectangular cartesian coordinates of the 

corresponding particle in the final state. Then the elastic deforma- 
tion is represented by particle-to-particle transformations 

(8-1) aF). 

Hence by the ordinary differential calculus 

(8 • 2) dx^ = /} da\ 

where 

fi _ aF) 

The classical theory of elastic bodies assumes that the deformations 8-1 
are “infinitesimal.’^ Such crude approximations have been found 
inadequate in some investigations on thin plates and shells. t As a 
result the finite deformation theory is beginning to be used in engineer- 
ing problems. In what we shall have to say we shall make the restric- 
tive assumptions of the classical infinitesimal theory only toward the 
end of the chapter. 

Let A. and X be defined as the column matrices of three elements i 
A = II II and Z = || a;* || . 

Similarly the differential matrices are dA = || da® || , dX = || dx^ || . 
Define the square matrix F by 

t See the Kdrman Anniversary Volume, California Institute of Technology, 1941, 
for various, papers and other references. 
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i.e., the matrix of the partial derivatives of the deformation 8 • 1. Then 
the differential relations 8 • 2 can be written in matric form as 


(8-3) 


dX dA, 


Definition. The adjoint M* of a matrix M is the matrix obtained 
from M by interchanging the rows and columns of M. 

Thus .4* and AT* are row matrices while 


jP* = 



do;® 

da;® 

dal' 

dal’ 

dal 


dx2^ 

da;® 

/O 

da^’ 

da® 



da;® 

da® 

da®’ 

da® 


It can be proved by a routine procedure that (MiAf 2 )* = M^M\. 
In other words, the adjoint of the product of two matrices is the 
product of their adjoints in the reverse order. For example, 

dX^ = dA'*F*. 

Hence the square of the differential line element in the final state of 
the medium will be 

(8-4) ds\ = dA^F^FdA 

since 

ds%: = {dx^y + {dx^y^ H- {dx^y = dX* dX. 

By a direct computation it can be shown that the matrix 

= II ypa 




(8-5) 
where 

^=1 c>a* da' 

Note that = ypji. This is expressed by saying that F*F is a symmetric 
matrix. 

The square of the differential line element in the initial state is 
(8-7) dsl = dA'^ d A (= dA'^I dA, where I is the unit matrix) 

and hence with the aid of 8-4 we find the formula 
( 8 - 8 ) dsi - dsl = dAr^F^'^F - I) dA. 

On defining the matrix, called the deformation or strain matrix j 

(8-9) H = 

we find 

(8-10) dsS: - dsl =^2 dA’^H dA. 
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Now, if dsx — dsA for all particles of the initial state and for all 
dA, we have, by definition, a rigid displacement of the medium from 
the initial state to the final state. A glance at 8*10 shows that a 
necessary and sufficient condition that the change of the medium from the 
initial to the final state he a rigid displacement is that the strain matrix 
H he a zero matrix. In other words, when H = O, the medium is not 
deformed or strained but is merely transported to a different po- 
sition by a rigid displacement. This property then justifies the 
terminology “strain matrix H” since H measures, in a sense, the 
amoimt of strain or deformation undergone by the medium. It is clear 
from definition 8-9 that the strain matrix H is a symmetric matrix. 
Let nij be the elements (more commonly called components in elasticity 
theory) of the strain matrix Hj i.e., -S’ = || || . On using result 

8 • 5 and definition 8 • 9, we see that 



Initial state Final state 

Fia. 8*1. 


Let = X* a^, then 

ba;* Zyu* 

== -}- S - 

baf ^ 

Define the matrix 

bw* 
ba* 

and obtain the relation 

F = F + /. 

Hence, from definition 8-9, for the strain matrix H we obtain 

H ~ il(V* 4 - r)(.V + /)_/]. 
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Expanding the right-hand side gives 

In the classical infinitesimal theory of elasticity only first-degree terms 

in — are kept. Hence, to the degree of approximtion cmtempMed iy 
00 ’ 

injinMnd Ikory of elastic deformtm, h strain rndrix is given ly 

E^\{y*^y) 


in rdangukr cartesian coordinates. In other words, the compents 
ilij of H are given by the familiar 





t 


From the symmetry of i)<y there are thus in general five distinct com- 
ponents of the strain matrix in three dimensions while there are three 
components for plane elastic problems. 



PART II. TENSOR CALCULUS AND 
ITS APPLICATIONS 

CHAPTER 9 

SPACE LINE ELEMENT IN CURVILINEAR 

COORDINATES 


Introductory Remarks. 

The vague beginnings of the tensor calculus, or absolute differential 
calculus £is it is sometimes called, can be traced back more than a 
century to Gauss’s researches on curved surfaces. The systematic 
investigation of tensor calculi by a considerable number of mathema- 
ticians has taken place since 1920. With few exceptions, the applica- 
tions of tensor calculus were confined to the general theory of relativity. 
The result was an undue emphasis on the tensor calculus of curved 
spaces as distinguished from the tensor calculus of Euclidean spaces. 
The subjects of elasticity^ and hydrodynamics,^ as studied and used 
by aeronautical engineers, are developed and have their being in plane 
and solid Euclidean space. It is for this reason that wo shall be pri- 
marily concerned with Euclidean tensor calculus in this book. We 
shall, however, devote two chapters to curved tensor calculus in con- 
nection with the fundamentals of classical mechanics and fluid 
mechanics. 

It is worthy of notice that the tensor calculus is a generalization of 
the widely studied differential calculus of freshman and sophomore 
fame. In fact, as we shall see, a detailed study of the classical dif- 
ferential calculus along a certain direction demands the introduction 
of the tensor calculus. 

Notation and Summation Convention, 

Before we begin the study of tensor calculus, we must embark on 
some formal preliminaries including some matters of notation. 

Consider a linear function in the n real variables x, y, z, • • • , it? 

(9*1) ouaJ + iSy -h -b • * • -h \w. 

Define 

Oi = a, 02 = |S, as = T, • • • , On = X, 

= y, 3^ = Zj " = w. 

42 
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We emphasize once for all that x'-, x”- are n independent vari- 

ables and not the first n powers of one variable x. In terms of the 
notations of 9 • 2 we can rewrite 9 • 1 in the form 

(9 • 3) aix^ + a2X^ + azx^ + • • • 4- OnX”, 

or as 

(9-4) ^aiX*. 

The set of n integer values 1 to n is called the range of the index i in 
9*4. A lower index i as in ai will be called a subscript, and an upper 
index i as in x^ will be called a superscript. Throughout our work we 
shall adopt the following useful summation convention: 

The repetition of an index in a term once as a subscript and once as a 
superscript will denote a summation with respect to that index over its 
range. An index that is not summed out will be called a free index. 

In accordance with this convention then we shall write the sum 
9 • 4 simply as 

(9 • 5) UiX ^ . 

A summation index as z in 9 • 5 is called a dummy or an umbral, since it 
is immaterial what symbol is used for the index. For example, a^x^ 
is the same sum as 9 • 5. All this is analogous to the (umbral) variable 
of integration x in an integral 

J f(x) dx. 

a 

Any other letter, say y, could be used in the place of x. Thus 

r /(?/)% = f fix') dx. 

•J CL ^ CL 


Aside from compactness, the subscript and superscript notation 
together with the summation convention has advantages that will 
become evident later. 

As a further illustration of the summation convention, consider the 
square of the line element 

(9 * 6) d^ — dx? -I- dy"^ 4- d^ 

in a three-dimensional Euclidean space with rectangular cartesian 
coordinates x, y, and z. Define 

(9-7) = X, y^ y, y^ = z 

and 

511 = ^22 = 8zz — 1, 

512 = 521 = 5i 3 = Ssi = 523 = 532 “ 0. 


(9-8) 
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Then 9 • 6 can be rewritten 

(9-9) 

or again 

(9 • 10) ds^ = 8ij dy^ dy^ 

with the understanding that the range of the indices t and i is 1 to 3. 
Note that there are two summations in 9* 10 one over the index i and 
one over the index j. 

Let f(x^, x^, • • ■ , x^) be a function of n numerical variables z\ x\ 
• • - , x^; then its differential can be written 


d/ = ^ (fa* 
da;* 


with the UTidersto/ndiug that the swnwiciticm coTiveution has Peen extended 
so as to apply to repeated ^per scripts in differentiation formulas. We 
shall adhere to this extension of the summation convention. 

It is worth while at this early stage to give an example of a tensor 
and show the fundamental nature of such a concept even for ele- 
mentary portions of the usual differential and integral calculus. This 
will dispel, I hope, any illusions common among educated laymen 
that the tensor calculus is a very ‘‘highbrow^' and esoteric subject 
and that its main applications are to the physical speculations of 
relativistic cosmology. 


Euclidean Metric Tensor.^ ' 

In the following example, free as well as umbral indices will have 
the range 1 to 3 as we shall deal with a three-dimensional Euclidean 
space. Let 

(9 • 11) = ff(y\ 2/*) 

be a transformation of coordinates from the rectangular cartesian 
coordinates y^, y^, y^ to some general coordinates x^, x^, x^ not neces- 
sarily rectangular cartesian coordinates; for example, they may be 
spherical coordinates. The inverse transformation of coordinates to 
9*11 is the transformation of coordinates that takes one from the 
coordinates x^, x^, to the rectangular cartesian coordinates y^, y^, y^. 
Let 

(9 • 12) yi = g^x^j x^) 

be the inverse transformation of coordinates to 9-11. 
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Example 

Let y^, y^, 2 /® be rectangular cartesian coordinates and x^, x^, 
polar spherical coordinates. The transformation of coordinates from 
rectangular cartesian to polar spheri- . 

cal coordinates is clearly ^ | 


xi = V{y^Y + (y2)2 4. 


x^ =» COS“^ 


+ W + Wj 


= tan~^ f ^ 


■© 


a!® (a;^, 


\ / 

V 

Fig. 9 - 1 . 


The inverse transformation of coordinates is given by 

y^ = x^ sin cos x^ 
y^ = x^ sin x^ sin a? 
y^ = x^ cos x^. 

The differentials of the transformation functions in 9*12 may be 
written 


(9-13) 


Zyu* 


On using 9 • 13 we obtain, after an evident rearrangement, the formula 
(9 • 14) <is“ - Z ^ ^ dx“ d^. 

OX ox'^ 


If we define the functions Qapix^, x^, x^) of the three independent vari- 
ables x^, x^, x^ by 


(9-15) 


g^{x\ - g 


we see that the square of the line element in the general x^j 0 ^, # coordi- 
nates takes the form 

(9 • 16) ds^ = dx“ dxP 

This is a homogeneous quadratic polynomial, called quadratic differen- 
tial form in the three independent variables dx^, dx^, da^. 

Caution: Once an index has been used in one summation of a series 
of repeated summations, it cannot be used again in another summation 
of the same series. For example, ga^ dx^ dx^ has a meaning and is equal 
to gn(dx^y + g 2 z(dx^y -f gz&{de^Y, but that is not what one gets by 
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carrying out the double repeated summation in 9 ‘16. Expanded in 
extenso, 9-16 stands for 

CQ.^7^ “ gxi{dx^y + 2^12 dx^ dx^ + 2gri3 dx^ da^ 

^ ^ I + g^idx^y + 2gr23 dx^ dx^ + gzz{dx^y. 


The factor 2 in three of the terms in 9 • 17 comes from combining terms 
due to the fact that is symmetric in a and /S; a glance at the defi- 
nition 9-15 shows that 


and hence 


Octp = g^cL for each ot and /3 


gi2 — 9 ^ 21 , 9^13 = 9 ^ 31 , 9^23 = 9 ^ 32 - 

Now let x^f x^j ^ be any chosen general coordinates, not necessarily 
distinct from the general coordinates Let 

(9 • 18) a;* = F^{x^, x^, 

be the transformation of coordinates from the general coordinates x^, 
x^, ^ to the general coordinates x^, x^, x^. Clearly the differentials dx'°‘ 
have the form 

(9 • 19) dx^= — dxP^. 

ox^ 


Define the functions gysi^^j of the three variables x^, x^y x^ by 

(9*20) gys{^^, x^) = 9'c.^(x\ x^y 

Then, if we use 9*19 in 9 • 16, we obtain, with the aid of the de fini tion 
9 • 20, the formula 

(9-21) ^3^ f^x"^ (^x^^ 

which gives the square of the line element in the OCi j y coordinates. 
We have thus arrived at the following result: 

If x^y x^j and x'^, x^, x^ are two arbitrarily chosen sets of general 
coordinates, and if the transformation of coordinates 9 ’IS from the x^s 
to the x^s has suitable differentiability properties, then the coefficients 
gasix^i x^y x^) of the square of the line element 9-16 in the x^, x^, co- 
ordinates are related to the coefficients gap(x^, x^, x'^) of the square of the 
line element 9 '21 in the coordinates x^, x^, by means of the law of tj'ans- 
formation 9 • 20. 

In each coordinate system with coordinates x^, x^, x^, we have a 
set of functions gasCx^, x^, x^'), called the components of the Euclidean 
metric tensor (field), and the components of the Euclidean metric 
tensor in any two coordinate systems with coordinates ^ y y and 
x^y x^y x^ respectively are related by means of the characteristic rule 
9-20. 
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An analogous discussion can obviously be given for the line element 
and Euclidean metric tensor of the plane (a two-dimensional Eu- 
clidean space). We now have two coordinates instead of three so 
that the range of the indices is 1 to 2. For example, the line elements 
ds in rectangular coordinates ( 2 /S y^) 

d^ = {dy^Y -1- {dy^Y 

will become 

daY 

in general coordinates (x^, and the components of the plane Eu- 
clidean metric tensor Qa^ix^j x^) will undergo the transformation 

g^,(x\ 3 ?) = 5^ 5^' 


Exercises 


1. Find the components of the plane Euclidean metric tensor in polar coordi- 
nates x^') and the corresponding expression for the line element. 

W - V(!/‘)= + W, 


« x^ cos x^, j / 9 

^ ’ and < „ . / 


) 2/2 = sin 
Hence 


x^ 


\ V(y')“ + wv 


y 


QijXX ,x)- 




y 




X 


y 


Fig. 9*2. 


Therefore 

a;2) = cos^ x^ + sin^ =■ 1, 

gizix^, x^) = (cos a:2)( -x' sin x^) + (sin x^Xx^ cos = 0 = g-nix^, a;*), 
giXx\ x^) = ( -x^ sin a;2)2 + (a;i cos x-y = ix^Y- 

The line element ds^ = gap da^ dx^ in polar coordinates (a;^, a;®) is then 

ds^ = (da;i)2 + {x^)Kdx^)\ 

which in the usual notation is written 

= dr^ + r2 dff^. 


2. Find the components of the (space) Euclidean metric tensor and the expres- 
sion for the line element in polar spherical coordinates. 

Answer. 

gii{x^, a:2, a;^) = 1, gi^fx'^, x^, a^) = (a;OS gzXx^, x\ x^) = (xO^sin x^y 
and all other gi^x'-y x^, x^) = 0, so that 

cfe2 » (da;i)2 + (xi)2(dx2)2 + (x0'*(sin x^yidx^y. 



CHAPTER 10 


VECTOR FIELDS, TENSOR FIELDS, AND EUCLIDEAN 

CHRISTOFFEL SYMBOLS 


The Strain Tensor. 

Another interesting example of a tensor is to be found in elasticity. 
Let cO-j a?, be the curvilinear coordinates of a representative particle 
in an elastic medium, and let be the coordinates of the repre- 

sentative particle after an elastic deformation of the medium. 



Medium Deformed medium 

Fig. 10'1- 


Let 

' 1) dsj = da“ daP 

be the square of the line element in the medium, and let 

’ 2) dx'°‘ dx^ 

be the corresponding square of the line element in the deformed medium 
induced by the elastic deformation whose equations are 

x^ a?, a?). 

In terms of the coordinates x^, 3 ? the line element 10-1 can be 
written 


(10 • 4 ) 
where 


(10 • 5 ) 

On subtracting 
( 10 - 6 ) 


, _ da* 

corresponding sides of 10-2 and 10*4 we fin d 


— dsl = 2e«j3 dxP, 
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if we define by 


(10-7) - K,). 

Now are functions of the coordinates x^, re®. If then we calcu- 
late ds^ — dso in any other curvilinear coordinates x^j we would 

obtain by the method of the preceding chapter 


where 

(10-8) 


•— d^ = 2eap dx"* d^. 


Because of the characteristic law 10-8, €„p are the components of a 
tensor (field), and because €ocp in 10*6 is a measure of the strain of the 
elastic medium, eap are the components of a strain tensor. We shall 
have a good deal to say about the strain tensor in some of the later 
chapters. 


Scalars, Contravariant Vectors, and Covariant Vectors. 

We shall now begin the subject of the tensor calculus by defining 
the simplest types of tensors. An object is called a scalar (afield) if in 
each coordinate system there corresponds a function, called a com- 
ponent, such that the relationship between the components in (x^y x^j a;®) 
coordinates and (x\ x^, 0) coordinates respectively is 

(10-9) x^y x^) = .s(x^, x^y 0). 

An object is called a contravariant vector field (an equivalent termi- 
nology is contravariant tensor jield of rank one) if in each coordinate 
system there corresponds a set of three functions, called components, 
such that the relationship between the components in any two coordi- 
nate systems is given by the characteristic law 

— 

(10*10) $*(^S x^j ^^) = ^ix^i xP-y a:®) r— r* 

ox 

An object is called a covariant vector field (an equivalent terminology 
is covariant tensor field of rank one) if in each coordinate system there 
corresponds a set of three functions, called components, such that the 
relationship between the components in any two coordinate systems 
is given by the characteristic law 

( 10 * 11 ) x^, 0) = Voeix^) x^y x^) — — 1 * 

O0 

It is to be noticed at this point that the laws of transformation 
10*10 and 10*11 are in general distinct so that there is a difference^ 
between the notions of contravariant vector field and covariant vector 
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field. However, if only rectangular cartesian coordinates are con- 
sidered, this distinction disappears. It is for this reason that the 
notions of contravariant as distinguished from covariant vector fields 
are not introduced in elementary vector analysis. That the char- 
acteristic laws 10*10 and 10*11 are identical in rectangular cartesian 
coordinates follows from some calculations leading to the result 


( 10 * 12 ) 


c)rc“ bx^ 


between rectangular coordinates a;®) and any other rectangular 

coordinates x% x®). 


Tensor Fields of Rank Two. 

In all three objects — scalar fields, contravariant vector fields, and 
covariant vector fields — there are components in any two coordinate 
systems, and the components in any two coordinate systems are re- 
lated by characteristic transformation laws. We have to consider 
other objects, called tensor fields (of various sorts) , whose components 
in any two coordinate systems are related l)y a characteristic trans- 
formation law. To shorten the statements of the following definitions 
we shall merely give the characteristic transformation law of com- 
ponents. 

Covariant Tensor Field of Rank -Two. 

— ^x^ 

(10*13) W(^S == ^'0 5^5^' 


Contravariant Tensor Field of Rank Two. 

br“ 

(10*14) i“^(x\ x2, X®) = x“, X®) 


Mixed Tensor Field of Rank Two. 


(10*15) 


t% (xS X®, X®) = ii(xb X®, X®) 


c)x“ bx^ 
bx^bx^ 


Again because of relation 10*12, the difference between the above 
three types of tensor fields is non-existent as long as one considers 
ordy rectangular cartesian coordinates. 

It is worthy of notice at this point that tlui indices in the various 
characteristic transformation laws tell a story which depends on 
whether the index is a superscript, called contravariant index, or a 
subscript, called covariant index. 

Before proceeding any further with the dev(dopinent of our subject, 
it would be illuminating to have some examples of vector fields and 
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other tensor fields. Perhaps the most important example of a contra- 
variant vector field is a velocity field. Suppose that the motion of a 
particle is governed by the differential equations 

dx^ 


where t is the time variable. If x^ = a;®) is a transformation 

of coordinates to new coordinates then 


dx^ 

dt 


'bx'^dx^ 
bx^' dt 


— x‘). 


Thus the components x^, x^) of the velocity field in the co- 
ordinates are related to the components x^, x^) in the x^ coordi- 

nates by the rule 



^%x^, x^, x^) = — ^^‘{x\ x\ x^), 


which is precisely the contravariant vector field rule. 

If s(x^, x^, x^) is a scalar field, then the “gradient’^ - — ; are the com- 

ox"^ 

ponents of a co variant vector field. An important example of a scalar 
field is the potential energy of a moving particle. 

We gave two examples of a covariant tensor field of rank two: the 
Euclidean metric tensor and the strain tensor. A little later in the 
chapter we shall give an example of a contravariant tensor field of rank 
two, the associated with the Euclidean metric tensor 

As an example of a mixed tensor field of rank two, we have the 
mixed tensor field with constant components 

dp = 0 if <x 9^ 13 
= 1 if a = jS 

in the x* coordinates. But 


Hence 


ix\ x^ 


bx^ bx^ 
bx^ bx^ 


50 (x^j x^, x^) =0 if <x 9 ^ ^ 

= 1 if O' = /3. 


In other words, not only are the components constant throughout space, 
but they are also the same constants in all coordinates. 

One of the first fundamental problems in the tensor calculus is to 
extend the notion of i^artial derivative to the notion of covariant deriva- 
tive in such a manner that the covariant derivative of a tensor field is 
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also some tensor field. It is true that, if one restricts his work only 
to cartesian coordinates (oblique axes), then the partial derivatives 
of any tensor field behave like the components of a tensor field under 
a transformation of cartesian coordinates to cartesian coordinates. 
For example, suppose that (x^, x^) are cartesian coordinates and 

(x^, x^y x^) are any other cartesian coordinates; then it can be shown 
that for suitable constants a} and a* 

(10 ■ 16) 0 = ajx^' + 

is the transformation of coordinates taking the cartesian coordinates 
X* to the cartesian coordinates x^. Hence 


(10-17) 

a set of constants, and so 
(10-18) 


da;*’ 



da;* da;* 


= 0 . 


Now, if ^^x^j x^y a;®) are the components of a contra variant tensor 
field, then 


(10-19) 


x^y 0) — ^“(a;S x^y a;®) 


d:r^ 

da;“ 


On differentiating corresponding sides of 10 - 19, we obtain 


( 10 - 20 ) 


d|^ d^ da;'* d;c^ ^ d^a;^ da;'* 

di^’ “ d^^ ^ ^ da;“ da;^ to ’ 


But both x^ and x^ are cartesian* coordinates. Hence on using 10 - 18 
in 10 ■ 20 we obtain 

d|^’ d^ bx^ dx* 
dx'* dx* dx“ ‘ 


( 10 - 21 ) 


dx* 


2>r 


which states that the partial derivatives — behave as though they were the 

ox^ 

components of a mixed tensor field of rank two and this under a trans- 
formation from cartesian coordinates to cartesian coordinates. 

The presence of the second derivative terms in 10 • 20 in curvilinear 

d|" 

coordinates x* shows that the are not really the components of a 

tensor field. So the fundamental question arises whether it is possible 
to add corrective terms C% (all zero in cartesian coordinates) to the 

partial derivatives so as to make 


dr 

dx^ 




( 10 - 22 ) 
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the components of a mixed tensor field of rank two for all contravariant 
fields The answer to this question is in the aflfirmative, and the 
possibility of the corrective terms depends on the existence of the 
Euclidean Christoffel symbols.® 

Euclidean Christoffel Symbols. 

We saw in the previous chapter that the element of arc length 
squared in general coordinates takes the form 

(10*23) ds® = rc®, x^) dx^ dx^, 

where, in the present terminology, are the components of a covari- 
ant tensor field of rank two, called the Euclidean metric tensor. Now 


(10*24) 

Define 

(10*25) 



9ii, 

912, 

9 13 

9 = 

^ 21 , 

922, 

923 


9sh 

932, 

933 


- 


9^ 0. 


Cofactor of in g 

■ " ' ' ■ P.I— ■ t 

9 


As the notation indicates, it can be proved that the functions g“^ are the 
components of a contravariant tensor field of rank two with the following 
properties: 

(10-26) 

(equals 0 if o: jS and 1 if a = jS). 

Define the Euclidean Christoffel symbols a:®, x^) as follows: 

— -I- cur _ 


(10*27) 


rLj 3 (a:S a;®, x^) = ^g 


Vda;“ ■ hx^ / 

Since the law of transformation of the components and are 
known, one can calculate the law of transformation of the r«^(ajS x®, a^) 
under a general transformation of coordinates x^ = ffipd-, ^c®, a;®). 

Let Qctff and g*^ be the components of the Euclidean metric tensor 
and its associated contravariant tensor respectively in the x^ coordi- 
nates. Then if we define 

(10 ■ 28) ^ 

we can prove by a long but straightforward calculation that the 
rL /3 (^^ 2 :®) are related to the r^^(:r^, a;®, 5®) by the following famous 

transformation law:^ 

— - ^x^ drc" d®x^ bx* 
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In the previous chapter we saw that 
(10.30) = 

where the y^s are rectangular cartesian coordinates. Consequently, 
if the a;’s are cartesian coordinates, it follows that all the components 


9txp(^^} 3*^6 constants. In other words. 




= 0 in cartesian 


coordinates x^. We have then immediately the important result that 
the EztcUdean Christoffel symbols x^, x^) are identically zero in 

cartesian coordinates. 

If the y^ are cartesian coordinates and the are general coordinates, 
one can calculate the Euclidean Christoffel symbols Tl,p{x^, x^, a:®) 
directly in terms of the derivatives of the transformation functions in 
the transformation of coordinates 

= fKy^, 2 /^, 2 /®) 

and in the inverse transformation of coordinates 


2 /* = x^, x^). 

Since all the Euclidean Christoffel symbols are zero when they are 
evaluated in cartesian coordinates 2/*, it follows immediately from the 
transformation law 10*29 that the Euclidean Christoj^el symbols in 
general coordinates x^ are given by the simple formula 

b^y^ Z>x^ 


(10-31) 


33®, X^) = 


da;“ Z>t/ 


This formula is often found to be more convenient in computations 
than in the defining formula 10 • 27. 

Caution: The Christoffel symbols are not the components of a tensor 
field so that z, a, and /3 are not tensor indices; i.e., i is not a contra- 
variant index and a., jS are not covariant indices. 

The concepts of tensor fields and Euclidean Christoffel symbols can, 
by the obvious changes, be studied in plane geometry — two-dimen- 
sional Euclidean space. Since we have two coordinates for a point in 
the plane, all components of tensors and the Christoffel symbols will 
depend on two variables, and the range of the indices will be from 
1 to 2 instead of 1 to 3. Thus the Euclidean Christoffel symbols for the 
plane uhll be 

(10-32) r*,(a:S x‘) = ^g^(x', ^ 


in terms of the Euclidean metric tensor g^pCx^, x®) for the plane. The 
alternative expression in terms of the derivatives of the transformation 




Exercise 


Compute from the definition 10 -32 the Euclidean Christoffel symbols for the 
plane in polar coordinates. Then check the results by computing them from 10 -SS. 

1 

HU: Use results of exercise 1 of Chapter 9 and find = 1, f * 0, f • 

(IJ 

inm. 




1 nl 


> Fli = 0, r|i = 0, Fij o Fji s 0, 122 = 0. 
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Covariant Differentiation of Vector Fields. 

Having shown the existence of the Euclidean Christoffel symbols, 
we are now in a position to give a complete answer to the fundamental 
question — enunciated in the previous chapter — on the extension of 
the notion of partial differentiation. We shall now prove that the 
functions 

(11 • 1) ' ■ - + nA a^) 


4- x\ (c®) a;®) 


are the components of a mixed tensor field of rank two, called the co- 
variant derivative of This result holds for every differentiable con- 
travariant vector field with the understanding that the x®, a;®) 

are the Euclidean Christoff el symbols. We shall use the notation P,t 
for the covariant derivative of By hypothesis we have 


(11.2) 


x^f a;®) = a:®, a:®) 


If we then differentiate corresponding sides of equation 11.2 we 
obtain, by the well-known rules of partial differentiation, 

bt ^x'" ^ b ’a:* b.r'" 

^ dr “ baf d.x“ d? ^ d^dr di“ ' 


We also have 
(11-4) 


da:* da:^ do:’ 


d®a-^ bx^ 


^j.y 


If we multiply corresponding sides of 11-4 by and sum on j we obtain 

(1 1 ■ 5) cvr dx’' ^ dx" dx> di“ ^ 


f -I- 

A ja? A 5- _QJ X .-X “ ? 


on using 


r = 


-.da:* 


in the first set of terms and 


50 
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in the second set of terms. Since X and ix are summation indices, we 
can interchange X and fx in the second derivative terms of 11-5. 
On carrying out the renaming of these umbral indices, we can add 
corresponding sides of 11.3 and 11.6 and obtain 

SlJL< 


( 11 - 6 ) 

Now 

(11.7) 

where 


'2)^® - 
I p 


U' = (; 


I •pX 

5^ + 



bx 


<>^x‘^ Z>x^‘ b:r*\ 

^ \bx^ ba^ b;r“ ba;^ b:r^' b:r“ 


Z>x^ V 

= 51, 

ba:'* b^« 


i _ / 1 if i 
“ (0 if i 


= a, 

9^ CL. 
X 


On differentiating 11*7 with respect to a; , we obtain 


( 11 * 8 ) 


b^x* ba;^ b:r* b^a;^ b:r^ 


bx^ bx^ bx“ ba/* bx^ bx“ bx^ 
and hence 11-6 reduces to 


= 0 


(11.9) 


-i - 

_±_ + r» - 
bx“ 


li' = (■ 


TlX 


But 11.9 states that the functions 

'bt 


bx^ 


+ rsu? 


are the components of a mixed tensor field of rank two. This completes 
the proof of the result stated at the beginning of this chapter. 

By a slight variation of the above method of proof, it can be es- 
tablished that the functions 


are the components of a covariant tensor field of rank two whenever are 
the components of a covariant vector field, called the covariant derivative of 
^i. As before, are the Euclidean Christoffel symbols. We shall use 
the notation for the covariant derivative of $». 


Tensor Fields of Rank r = p q, Contravariant of Rank p and Co- 
variant of Rank q. 

It is convenient at this point to give the definition of a general tensor 
field. As in the case of tensor fields of rank two, the definition will be 
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clear if we give the law of transformation of its components under a 
transformation of coordinates. 






( 11 . 11 ) Tt:::r.(ix\ xs x^ x^)^^ 




bx"*^ bx**^ 

^ b^ ‘ ■ ‘ 


We are now in a position to consider some problems that arise in 
taking successive covariant derivatives of tensor fields. However, we 
must first say a word or two about the formula for the covariant deriva- 
tive of a tensor field. If ‘ ' are the com'ponents of a tensor field of 

rank p + q, contravariant of rank p and covariant of rank g, then the 
functions defined by 


( 11 . 12 ) 


{ 


rpoci 


CCp 

Pq* y 


^nrrai • • • ocp 

pi* * ^ pq . -pofi rp<rot2 
/5i- 


•arp 

Pa 


+ 


+ xioip rpoti 
<ry-^ pi 


ccp-i<r -po- nfion • - 

Pq a‘p2 * 


QJp’’ 

^ Pa 


pier rpoti • • • ocp 
Pay-^ Pf • • Pg-ia 


are the componerds of a tensor field of rank p + q + 1, contravariant of 
rank p and covariant of rank q + 1. 

Tpl'. ^ will be called the covariant derivative of \ .* The above 
result ^ stating that the covariant derivative of a tensor field is indeed 
a tensor field can be proved by a long but quite straightforward calcu- 
lation analogous to the method of proof given for the co variant deriva- 
tive of vector fields. 

Since the covariant derivative of a tensor field is a tensor field, we 
can consider the covariant derivative of the latter tensor field, called 
the second covariant derivative of the original tensor field. In symbols, 
if is the original tensor field, we can consider its second suc- 

cessive covariant derivative 


rjioti • • • cLp 

/3i — fit, y. S' 

The fundamental question arises whether covariant differentiation is a 
commutative operation^ i.e., whether 

(11 ■ 13 ) Tff . : : z y, s = T%\: z. s. y • 


The answer is in the affirmative? In fact, since all the Euclidean Chris- 
toffel symbols are zero in cartesian coordinates, the partial derivatives 
of all orders of the Christoffel symbols are also zero in cartesian coordi- 
nates. Hence if *TZ'.y.%l{y\ 2/^ 2 /®) are the components of T^'.WZ in 
the cartesian coordinates y^, we find with the obvious repeated use of 
formula 11.12 that 




OCp 

pQ* T. S 


(.y\ y^) = 


>.2^ moci ' • • ctp 

by'^ by^ 
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In other words, we have the result that successive covariant derivatives 
of a tensor field reduce to partial derivatives of the tensor field whenever 
the tensor field and the operations are evaluated in cartesian coordinates. 


Properties of Tensor Fields. 

Perhaps the most important property of tensor fields is the following: 

If the components of a tensor field vanish identically (or at one point, or 
at a set of points) in one coordinate system, they vanish likewise in all 
coordinate systems. This result follows immediately on inspecting the 
law of transformation 11-11 of the components of a tensor field. 

If, then, one can demonstrate that a tensor equation 




holds good in one coordinate system, it will necessarily hold good, 
without further calculation, in all coordinate systems. For example, 
consider the covariant derivative Py, ^ of the Euclidean metric tensor. 
Now, since 


9ij, k 


da;* 




and since the gij are constants when evaluated in cartesian coordinates, 
we have = 0 in cartesian coordinates, and hence the tensor equation 
Py, i = 0 holds in all coordinates throughout space. 


Exercises 

1 . Prove that the covariant derivatives of g'^ and of 5) are zero. See equation 
10-26 for the definition of the tensor The mixed tensor 8} = 1 if i j and 
8} = 0 if i 7^ j. 

2. If Tfi is any tensor field, then show that Ta is a scalar field. Similarly, if 
is a mixed tensor field of rank throe, show that 7’“^ covariant vector field.® 

3. If the tensor field Tfl is defined by T“ = nap in terms of the two tensor 
fields X"'* and prove that the following formula (reminiscent of the differentia- 
tion of a product in the ordinary differential calculus) holds: 

7 “ + X®*^ liap, y. 



CHAPTER 12 

LAPLACE EQUATION, WAVE EQUATION, AND POISSON 
EQUATION IN CURVILINEAR COORDINATES 


Some Further Coucepts and Remarks on the Tensor Calculus. 


We remarked in Chapter 10 that, if s(a;S is a scalar field, then 


bs 

bx 


is a covariant vector field. So, to complete the picture of covariant 

bs 


differentiation, we can call — ^ the covariant derivative of the scalar field 
s{x^, x\ a®). 

For some discussions it is convenient to extend the notion of a tensor 
field. By a relative tensor field of weight w, we shall mean an object 
with components whose transformation law differs from the transforma- 
tion law of a tensor field by the appearance of the functional determi- 
nant (Jacobian) to the tyth power as a factor on the right side of the 
equations. If it? « 0, we have the previous notion of a tensor field. 
For example: 




and 





bx^ 

bx^ 

bx^ 

bx^ 

bx® 

bx2 

bx* 

bx* 

bx^ 

bi*’ 

bx® 

bx® 

bx® 

bx® 

bx^ 

bx* 

bx® 

bx^ 

bx^ 

bx^ 

bx^ 

bx* 

bx® 

bx^ 

bx* 

bx* 

bx^ 

bx* 

bx® 

bx® 

bx® 

y 

bx® 




s{x^, x®, x^) 


w 


bx^ bx^ 




are the transformation laws for a relative scalar field of weight w and 
a relative contravariant vector field of weight w respectively. A rela- 
tive scalar field of weight one is called a scalar density , a terminology 
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suggested by the physical example of the density of a solid or fluid. 
In fact, the mass m is related to the density function p{x^j of the 

solid or fluid by 

m = J" S J' p(x^, x^y x^') dx^ dx^ dx^y 

where the triple integral is extended over the whole extent of the solid 
or fluid. 

Another important example of a scalar density is given hy y/g, 
where g is the determinant of the Euclidean metric tensor g^^. In fact, 

(12 * 1 ) ^ • 


Let g — I g^p | , the determinant of the gap- By a double use of the 
formula for the product of two determinants when applied to 12*1, 
one can prove that 

( 12 - 2 ) S-h 


( 12 - 2 ) 


where 


stands for the functional determinant of the partial 


derivatives- " — • 
ox 

On taking square roots in 12 • 2 we obviously get 


(12*3) 


/- /- 


which states that ■%/ gr is a scalar density. 

The g enters in an essential manner in the formula for the volume 
enclosed by a closed surface. In fact, the formula for the volume in 
general curvilinear coordinates x^ is given hy the triple integral 

(12*4) V = fff^ dx^ dx^ dxK 

This form for the volume can be calculated readily by the following 
steps. If 2/* are rectangular coordinates, then 

V = J'J'S dy^ dy^ dy\ 

Hence 

(12 -5) V = fSSJ dx^ dx^ dx^y 

where J stands for the functional determinant 

of the transformation of coordinates from the curvilinear coordinates 
x* to the rectangular coordinates y^. Now 

. , o nx ^ 2)2/» by^ 

x^ xfl) = g —3^, 
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and hence the determinant 

g = S'cjs 1 

is precisely equal to on using the rule for the multiplication of two 
determinants. In other words 

( 12 - 6 ) Vg = J, 

from which the formula 12 • 4 for the volume becomes clear. 

Since the formula for volume in general coordinates a;* is given by 
12-4 it follows that the mass w of a medium in general coordinates x* 
has the form 

m — x^j x^)\/g dx^ dx^ dx^, 

where po(x\ x^, x^) is an absolute scalar field that defines the (physical) 
density of the medium at each particle {x^, x"^, x®) of the medium. 
Clearly p(x^, x®) = po(x^, x^, 'x®)-\/gr is a scalar density and, since 

\/g = 1 in rectangular coordinates, has the same components as the 
density po{x^, x^, x®) of the medium in rectangular cartesian coordinates. 

Other concepts and properties of tensors will be discussed later in 
the book whenever they are needed. 


Laplace’s Equation. 

Let be the contravariant tensor field of rank two defined in 
Chapter 10, i.e., 

(12.7) _ Cofactor of in g 

9 

If \l/(x^, x^, X®) is a scalar field, then the second covariant derivative 
is a covariant tensor field of rank two. Now we can show that 
9"‘^^.a,fi a scalar field. In fact, 

?)rf“ 

( 12 - 8 ) = 


S^ = g^ 


(12-9) 


- dx*^ bx’’ 

’/'.a.P = ^,.,r ^ 


On multiplying corresponding sides of 12-8 and 12*9 and summing on 
OL and /3, we obtain 

-a /? , \u , 

g - g 

and hence the desired result 




on using the obvious relations 


(12-10) 


5x“ bx*^ 
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Let 


(12- 11) F(x\ x^, a?) = a:^ xF) 

for the arbitrarily chosen scalar field x^, x^'). W^e have just shown 
that F{x^, x^ f o;^) is also a scalar field. In rectangular cartesian coordi- 
nates, the Euclidean metric tensor has components equal to 0 for 
OL 9 ^ ^ and equal to 1 for ex. — Furthermore, we saw in Chapter 11 
that in cartesian coordinates, and hence in rectangular cartesian 
coordinates, ?/% 



(2/S 2/S 2/®) 


y^) 


The component of the scalar field F{x^j xp-, x^) in rectangular coordinates 
2 /* is then 


^^V(2/S 2/^, 2/^) 


or 

( 12 . 12 ) 


^^VCi/S y^, y^) y^) 2/^> 2/^) 

wy wy ^ wy 


the Laplacean of the function *^(y\ y\ y^). 

Hence the form of Laplace’s equation in curvilinear coordinates x* 
with the scalar field ^{x^, x^, x^) as unknown is given by 


(12. 13) 


x^) = 0, 


where g'°‘^{x^, x^, x^) is the contravariant tensor field 12.7 defined in terms 
of the Euclidean metric tensor g„^, and where x^) is the second 

covariant derivative of the scalar field ^{x^, x^, x^). 

If we write 12.13 explicitly in terms of the Euclidean Christofifel 
symbols r}A;(x^, x^, x®), we evidently have 


(12.14) g-'-CxS xS = 0 


as the form of Laplace’s equation f in curvilinear coordinates x*. 

It is worth while at this point to give an example of Laplace’s 
equation in curvilinear coordinates and at the same time review 
several concepts and formulas that were studied in previous chapters. 


t There is another form of Laplace's equation in curvilinear coordinates »* 
which is sometimes more useful in numerical calculations than 12- 14. It is given by 

1 

\7^ ^ proof see note 3 to Chapter 13. Similar remarks 

could be made for the wave equation and Poisson's equation since the Laplace 
differential expression occurs in them. 
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Let 2 /^, be rectangular cartesian coordinates and x^, x? polar 
spherical coordinates defined by the coordinate transformation 

iy^ — x^ sin x?- cos x?. 




(12 • 15) sin x?' sin x?^ 

\y^ = x^ cos x^. 

Clearly the inverse coordinate trans- 
formation is given by 

a;! = V (v')^ + W + (y^y 


X^ = COS" 


X? — tan 


(2/^ 

-(?)■ 


y + w + ( 2 /^)' 


Fig. 12*1. 


From the definition of the Euclidean metric tensor we find 
(12-16) g^ii = 1 , 6^22 = Qzz = {?yy (sin x^y^ and all other Qij — 0 , 

so that the line element ds is given by 
(12-17) ds^ = (dx^y + (x^y{dx^y 4 - (a;i)2(sin x^y{dx^y. 

Again, from the definition of the tensor we find 


(12-18) L = 7 -I 7 „ = 


, and all other — 0 . 


The Euclidean Christoffel symbols can now be computed in spherical 
polar coordinates; use either formula 10-27 or 10-31. They are as 
follows: 

(^22 = — rJa = — a:^(sin x^y, 

jl’ia “ rii = “> rfs = —sin x^ cos x^, 

(12-19) ( ® 

irfs = r|, =■ -, ria = ria - cot 
\and all other = 0 . 

On using 12-18 and 12-19 in 12-14 (r 22 , rfs are the only non- 
zero Christoffel symbols that are actually used), we find that Laplace’s 
equation in spherical polar coordinates x^, x^^ x^ is 

( by _ _I^ _dV . 1 dV 

)(bx^ 


(12-20) 


1- 1 dV 1 dV 

y {x^y (Px^y (a;^)2(sin x^y (Px^y 

2 cot x^ ^>\lf 

_l 1 - _| z_ ^ 0 

a;! bxi (x^y ba;2 


whenever the unknown function is a scalar field \p(x^, x^, x^). 
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Laplace’s Equation for Vector Fields. 

We now turn our attention to the related problem of considering 
vector fields whose individual components satisfy Laplace’s equation in 
rectangular coordinates. The question arises whether each component 
of the vector field will satisfy Laplace’s equation 12 - 14 in curvilinear 
coordinates x^. The answer is in the negative, as a little reflection 

will now show. To be specific, let the unknown be a contravariant 
vector field x^) with components 2/®) in rectangular 

cartesian coordinates. By hypothesis 


( 12 - 21 ) 


^ 2*^4 


+ 


+ 


= 0 . 




“/sf- 


Lbx' 




wy wy ' (^>y^y 

By practically the same type of argument used in deriving equations 
12-13, we find that the contravariant vector field x^, x^) in curvi- 
linear coordinates x^, x^, x^ will satisfy the system of three differential 
equations 

( 12 - 22 ) = 0 , 

where is the second covariant derivative of x^, x^). If we ex- 
pand 12-22 explicitly in terms of the Euclidean Christoff el symbols 
x^y x^) we find 

_ ^ _i_ p* ^ 

(12-23) -- 

, f crot 

a system of three differential equations in which all three unknowns 
occur in each differential equation. 

Wave Equation. 

The propagation of various disturbances in theory of elasticity, 
hydrodynamics, theory of sound, and electrodynamics is governed by 
the partial differential equation known as the wave equation. In 
rectangular cartesian coordinates y\ the wave equation is 

b^*u{y\ 2/2, 2/®, t) ^^*u{y\ y^, y^, t) 


f 511 

ba;“ 


nsKc - 


T^r 

<rr^ otP 


)r] - 0, 


(12-24) 


bi2 


i^y^y 

y\ y^y t) b^*^(2/S 2/^ y^, t) 

wy 


and hence in curvilinear coordinates x* 


(Py^y 
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where is the second covariant derivative of the scalar field 

u{x^j x^j x^, t). If we write 12-25 explicitly in terms of the Euclidean 
Christoffel symbols we obtain f 


(12-26) 


b‘^u(x^, x^, x^, t) 
_ 



bx^ bx^ 




It is to be observed that the right-hand side of 12 • 25, or equivalently 
of 12-26, is the Laplacean. If the x^ are spherical polar coordinates, 
Laplace’s equation has the form 12-20. Hence, immediately, we see 
that the wave equation has the following form in spherical polar co- 
ordinates 


(12 • 27) 


b'^u(x^j x^, x^, t) b^u 1 bhi 1 bhL 

bt^ (da;^)^ (bx^y (a:^)^(sin (bx^)^ 

2 bu cot x^ bu 
x^ bx^ (a^O"^ ^a;^ 


By exactly the same calculations as we used in obtaining Laplace’s 
equation for contravariant vector fields in curvilinear coordinates, we 
find that the wave equation in curvilinear coordinates x^ takes the follow- 
ing form whenever the unknown is a contravariant vector field ^*(a;^, x^j 
x^j t) that depends parametrically on the time t: 


(12-28) 

or in expanded form 


x^, x^, t) _ . 


bf^ 




(12-29) 


^^^Kx\ x\ x^ 0 ^ b^^^ _ r-" . r* 

bt^ ^ Lbx^^bx^ ^^bx^^ 


+ r; 






\ ba^ 


+ - r*.r, 


;«)r] 


on using the corresponding result 12-23 for Laplace’s equation. Note 
that 12-29 is a system of three differential equations for the three 
unknowns and and not just one differential equation satisfied 

by the three functions and 


Poisson’s Equation. 

As a final exercise in this chapter we consider Poisson’s differential 
equation. In rectangular cartesian coordinates y'^, Poisson’s equation is 

c)"V(2/h 2/^ y") ^y{y\y\y^) ^^*Hy\y^y^) 

(12-30) (by^y (by^y + (fyy^y 

= — 47r*<7(2/i, 2/^ 2/®). 

t For another form of the wave equation, see Exercise 1 at the end of the chapter. 
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On the right-hand side of 12 • 30, 2/®) is the component in rec- 

tangular coordinates 2/* of a scalar field 


(12 -SI) 


cr{x^, x^j x^) 


pix^f x^, x^) 

^/g 


•y 


where p{x^, x^, x^) is a scalar density. f Since -v / g is also a scalar density, 
c(x^, x^, x^) is obviously a scalar field. From the definition of gr as the 
determinant of the Qap, we see that in rectangular coordinates 


(12*32) 

Hence 


*g(.y\ 2 /S 2/®) = 


10 0 
0 10 
0 0 1 


- 1. 



and consequently 

(12*33) *<T{y\ y\ y^) = *p{y\ y\ y^). 


In other words, the scalar field a{x^j x^, x^) and the scalar density 
p(x\ x^, x^) have equal components in rectangular cartesian coordinates. 

On making use of our calculations for Laplace’s equation, we can 
derive corresponding results for Poisson’s differential equation. For 
example, Poisson's differe^itial equation in curvilinear coordinates x^ 
will be 


(12*34) 


g 





— 47rcr(a;', x®). 


whenever the unknown is a scalar field »/'(x^, x®, x®). As before, the are 
the Euclidean Christoff el symbols in the x^ coordinates. 


Exercises 


1. Show that the wave equation in curvilinear coordinates x*' with scalar uix\ 
x®, X®, t) as unknown can be written as 


c)®w(x*, X®, X®, 0 


1 

dx“ 


2. Show that Poisson’s equation in curvilinear coordinates x* with scalar \^(x^, x®, 
X®) as unknown can be written as 


1 



— 4Trcr(x^, X®, X®). 


t In most physical problems p(x^, x®, x®) = poCx*, x®, x^)y/g, where po(x^, x®, x®) 
is an absolute scalar field aud represents the physical density of a medium. 





CHAPTER 13 


SOME ELEMENTARY APPLICATIONS OF THE TENSOR 
CALCULUS TO HYDRODYNAMICS 


Navier-Stokes Differential Equations for the Motion of a Viscous 
Fluid. 

As an interesting application of the covariant derivative of a tensor 
field in hydrodynamics, we shall write the famous Navier-Stokes differ- 
ential equations in curvilinear coordinates. Let and y^ be rec- 

tangular cartesian coordinates, and let 

w* = u*{y^, y^ , 2 /^, t), the contravariant velocity com- 
ponents of a viscous fluid. 
t = time. 

P = P(2/S y^i y^f t)t pressure. 

P = p(y\ y^y y^, t)y density. 

At = coefficient of viscosity, a constant. 

II 

V - kinematic viscosity. 

X* = Z*(2/S y^y y^y contravariant vector components 
of body force per unit mass. 

Then the motion of a viscous fluid is governed by the four Navier- 
Stokes differential equations.^ 



(13-2) 


^u* _ / \ ^ bu® V b / bu^\ 


bi by® 



The last differential equation is the equation of continuity ^ which 
expresses the requirement that the mass of any portion of the liquid is 
conserved. For a non-viscous fluids /u = 0 and hence v = 0, the Navier- 
Stokes equations reduce to the Eulerian hydrodynamical equations.^ 

The expression within the parenthesis in 13.2 is the Laplacean of u'. 
If we then make use of the results of the previous chapter on the form 
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is a mixed tensor field of rank two. Hence, by the argument just com- 
pleted for 1 I , we see that the determinant \ 5| — 2€| | is itself a 
strain invariant. But, from the expansion 15 • 1, we see that 

(15-4) /a = iCl - 2e? | - 1 + 2/i + 8 / 3 ]. 

Formula 15*4 expresses I 2 as a linear combination of strain invariants 
ivith numerical multipliers. Hence obviously J 2 is itself a strain in- 
variant. 

On using the results 14*31 together with what we have just proved, 
we obtain the result 


(15*5) 'jrr~ ~ ^ ^-^2 — 8 / 3 , 

dV 

which gives the ratio of the element of the volume of a set of particles in the 
unstrained medium to the element of volume of the corresponding particles 
in the strained medium in terms of the three strain invariants /i, I 2 , and I 3 , 


Homogeneous and Isotropic Strains. 

Let us now discuss the mathematical description of a homogeneous 
strain. 

Definition of Homogeneous Strain. A strain is homogeneous if 
the corresponding strain tensor has a zero covariant derivative, 
i.e., y = 0. 


In rectangular coordinates, the condition reduces to = 0 since all 

ox^ 

the Euclidean Christoffel symbols are identically zero in rectangular 
coordinates. In other words, the strain terisor components in rec- 
tangular coordinates are constants for a homogeneous strain. 

It readily follows from the definition of the strain invariants Ix, I2, 
and I 3 that for a homogeneous strain 




0 


in rectangular coordinates and hence in all coordinates. (Keep in mind 
that Iij I 2 , and Iz are three scalars and not the three components of a 


covariant vector.) Hence, for a homogeneous strain, is a numerical 

dV 

constant for all coordinates. We thus arrive at the important result 
that for a homogeneous strain 


(15*6) 


Vo 

V 


V l - 27i + 4/2 — 8/3 = a constant. 


This constant is the same for all unstrained^’ volumes Vq and their cor- 
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responding strained^ ^ volumes V, and is independent of the coordinate 
system. 

For the special case of a homogeneous strain which is also isotropic at 
each point, we shall have 

(15-7) = 

or what amounts to the same thing 

(15*8) (jg^, the Euclidean metric tensor.) 

Since the strain is homogeneous, we have .y = 0. We also have (see 
the end of Chapter 11) 

QuS, -y = 0 . 

Hence e in 15 • 7 and 15 Sis a constant scalar field, a numerical constant 
that is independent of position and the coordinate system. 

An example of an isotropic homogeneous strain is found in an iso- 
tropic medium subjected to uniform hydrostatic pressure, i.e., an 
isotropic medium subjected to the same pressure in all directions. 
Since 15 • 6 holds and 

1 _ 2/i + 4/2 - 8/3 - I - 26^ I , 

we see by an evident calculation that for an isotropic homogeneous strain 

(15-9) r “ 


The constant scalar e for an isotropic homogeneous strain is given 
by the formula 


(16*10) 



in terms of any one volume before and after deformation. 

In the usual approximate theory (usual theory of elasticity) higher 
powers of e than the first are neglected; see Chapter 14. Since 


we have 
(15*11) 
and 
(15*12) 


(1 — 2€)^ = 1 — 3e, approximately, 


Vo 

V 


= 1 — 36, approximately, 


IF- Fo 
* 3 F 


lAF lAF 


A Fu n d am ented Theorem on Homogeneous Strains. 

We shall now outline the proof of the following theorem. A necessary 
and sufficient condition that a strain he homogeneous is that, in terms of 
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unstrained cartesian coordinates “s and strained cartesian coordinates y*, 
the deformation is linear, i.e., 

(15-13) “2 = ^Ai 

Outline of Proof. 

From the definition of the strain tensor eap, we have 

hpq Gpg. ~~ 2€pg. 

Since gpq,r ~0, and since under our '‘necessity hypothesis” = 0, 
we obtain the vanishing of the co variant derivative of hpq (x). But 
by definition 

C^) “ (ai3®)( 

SO that the x* are the independent variables and the “a are the de- 
pendent variables. Expanding the covariant derivative in hpq^r (.x) — 0, 
and rearranging, we find 


(15-14) 

^,pr a^q ® a — 

where 


(15-15) 

« d“ap „ 

Since 


(15-16) 

“a,, = '‘a„, 


the right side of 15 • 14 must be symmetric in p and r. Hence 
(15-17) oe/gC “a,pr g = — ~ oc0p ^a^qv 


s _ cc fi 
a.r — a/jC »,r O.gp* 


On equating the corresponding sides of 15-14 and 15*17, rearranging, 
and interchanging p and q and jS and a, we obtain 

/IK - 1 0\ Sn „ «yTt ^<-1 _l_ «/ir Sn _ 1 _ _/» 


(15-18) cts^**a^pr ^a^q — — “a.p ^a^q + ar a^p + „pp a^r ««p- 
Adding corresponding sides of 15-17 and 15-18 there results 
(15 • 19) 2«^c “a pr ^a^q = — “o.r ^a^q “a p ^ “«.r 


Recalling that the are functions of the unstrained coordinates *a, we 
find 


1 ^«|8< 

2Le>"'a 


C b^.yC 


a^r ^Oi^q "^a^pt 


(15 • 20) a/3f' 
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Multiplying corresponding sides by and summing on q, we obtain 


(15*21) 


« _ if 

o,pr 2Lc>^a 



Finally, if we multiply corresponding sides of 15 • 21 by (f’^, we arrive 
readily at the interesting result 

(15*22) ®,pr ~ 'y«F(u) 

where Euclidean Christoff el symbols based on the Euclidean 

metric tensor cc^^^ia) in the unstrained coordinates ^a. 

Now from the definition 15*15 of “a,pr and from the vanishing of the 
Euclidean Christoffel symbols r}&(x) and }&r(o) when evaluated in 
“strained” cartesian coordinates and “unstrained” cartesian co- 
ordinates % respectively, we see that 15 • 22 reduces to 


(15*23) 




This implies that the deformation is linear, i.e., of type 15*13. 

To prove the converse part of the fundamental theorem on homo- 
geneous strains, we have by hypothesis that the deformation, or strain, 
is given by a linear transformation 15*13 in cartesian coordinates 
and y^. Now 

hpq(^x) = a!/3c(ci) ( a^p)^^a^q^ 


and are constants in cartesian coordinates, 
have 


d2/p 




p 


From 15*13 we 


and hence the components *hpq(y) in cartesian coordinates y* are 
given by 

*hpq(y') — A.p 
a set of constants. Hence 

d?/’- 

But this condition implies that the covariant derivative 

C'^) ~ n 

and hence the covariant derivative epq^r (^) = 0. In other words, the 
strain is homogeneous, and the proof of the theorem is complete. 


Variation of the Strain Tensor. 

In preparation for the subject matter of the next chapter, we need 
to consider deformations that depend on an accessory parameter 
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which in dynamical problems can be taken as the time t. So let the 
coordinates of a representative particle in the strained medium 
depend on the coordinates of the corresponding particle in the 
unstrained medium and on the accessory parameter t. Let 


(15-24) 


= 


U 


dt 


be the partial differential of a:Mn t. If f'. \ \ (x) is any tensor field in the 
strained medium, define bfW. {x) by 

(15 • 25) 8fv/. {x) = /: : 

where is the covariant derivative of f'.W. Clearly, if the are 

cartesian, 8f\ \ \ = T>f\ \ Moreover, 8f{x) = D/(a;) for a scalar f{x) in 
general coordinates x^. To have a well-rounded notation, define 
= Da:’’ and refer to bx'^ as the virtual displacement vector. If 
x^(}a, ^a, ^a, t) have continuous second derivatives, then from the 
commutativity of second derivatives 


Hence 




b(8x^) ^(5x0 


2)*a 




(a,3J**)* 


OX 


which implies the tensor equation 
(15-26) 5(dx*) = (5x*'),« dx^. 

Obviously dgra = 0, since grs,t = 0. Hence the above tensor equation 
may be written 

(15-27) 8(dxk) = (5 xa ;),« 

Since 8d(ja) — 0, an evident calculation using 15-26 shows that 
(15-28) 5(-a^) = -^a„(5x“).^. 

Recalling that 

hpq(x^ ~ o:/3<^(*^)( ®,p)(^^,a) 

and applying formula 15 • 28 we find 

bhpq = a0^\l fp^a^Q H- a^p ^a^j(fix ),<]r3l7 

which can be put in the convenient form 
(15-29) dhpq = -/iy(5x^),p - hp{8xf),q. 

From this, and from the definition of the strain tensor and the re- 
lated formula 

24 , 
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we arrive at ihe fundamental fomida for the variaim of Ote strain tmor. 
(15'30) iejjj ( 3 ) = IC(J3j),p + (53p)^J - [ 65 (S 3 ,)|P + ig 
This formula becomes 

(15-31) kg,[x) = ^[K).p + («3p),J 

within the ajiproximations of the usual approximate theory qf elasticity. 

Returning to our finite deformation theory, we define a rigid virtud 
disphemerti by the condition S{ds^) = 0 . On using formulas 15-26 
and 15-27 in an evident calculation, we find 

(15 - 32) S{d^) = [( 81 J ^ + (Sij) J dsf dxP 

for any virtual displacement, rigid or not. Hence the virtual displace- 
ment vector must satisfy Killing’s differential eguations for a rigid 
virtual displacement 

(15-33) (53«)j3 - 1 - {Sx^)^ = 0. 

For the sake of completeness, we shall write down the formula 
(without giving the derivation) for the variation of the Lagrangean 
strain tensor pjij under an arbitrary virtual displacement. 

M “ |[(to«),/s + (5*^) J pf a/. 


Exercise 

Calculate the three fundamental strain invariants for a homogeneous isotropic 
strain. Hint: since they are constants, calculate them in rectangular cartesian co- 
ordinates. 
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STRESS TENSOR, ELASTIC POTENTIAL, AND 
STRESS-STRAIN RELATIONS 


Stress Tensor. 

Let S be the bounding surface of a portion of the elastic medium 
in its strained position. The surface element of S may be described by 
means of the covariant vector dSr, 

(16 • 1) dSi = V? a?®), d/S 2 = -s/g d(a^, x^), dS^ * -s/g d(x\ x^), 


where 

(16-2) 


a;®) = 


bu bv 
bx^ bx^ 
bu bv 


du dv 


and Uf V are the surface parameters so that the parametric equations 
of the surface S are given by x^ = f^(u, v). In rectangular coordinates 
and in the usual notations x, z, the components of the covariant 
vector dSr are given by 

d/S» = d(2/, z), dSy “ d{Zj x), dSg = d(x, y). 


Exercise 

Prove that dSr is a covariant vector under transformations of the 
coordinates x^. Hint: use the fact that 'y/g is a scalar density. 

Let d/S be the magnitude of the surface element, i.e., 

(dS)^ = dS^ dSff, 

Before we introduce the notion of a stress tensor we must define a 
stress vector. A stress vector is a surface force that acts on the surfa>ce of 
a volume. An example of a surface force is the tension acting on any 
horizontal section of a steel rod suspended vertically. If one thinks of 
the rod as cut by a horizontal plane into two parts, then the action of 
the weight of the lower part of the rod is transmitted to the upper part 
across the surface of the cut. A hydrostatic pressure on the surface of 
a submerged solid body provides another example of a surface force. 

There are other kinds of forces called body, volume, or mass forces, 
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i.e., forces that act throughout the volume. As a typical example of a 
mass force one can take the force of gravity, pg AF, acting on the mass 
contained in the volume AF of the medium whose density is p, and 
where g is the gravitational acceleration. 

A stress tensor is defined implicitly hy the reloMon 

(16-3) F-dS^T^^dS^, 

where F^ is tlw stress vector acting on the surface element dSr. 

Let us now consider a virtual displacement of the strained medium 
corresponding to the accessory parameter t. The virtual work of the 
stresses across the boundary S is 

(16-4) J^y'F^ 8x^ dS = dSc, = dV, 

a volume integral extended over the volume F bounded by and ob- 
tained by Green's theorem or generalized Stokes' theorem in curvi- 
linear coordinates. 

If there are mass forces (M’’ per unit mass) acting on the medium, 
the virtual work of these mass forces is 

dV, 

V 

where p is the mass density. Hence, the virtual work of odX the forces 
acting on any portion of the medium is 

(16 • 5) fjrSXST^ + pM^)8x^ H- J dV, 

We shall now adopt the 

Physical Assumption op Equilibrium: The virtual work of all the 
forces acting on any portion of the medium is zero for any rigid virtual 
displacement. 

In particular, the translations, characterized by (8Xp)^q = 0, are 
rigid virtual displacements, and so we must have the condition 

(16 • 6) + pM^) 8x^ dV = 0. 

Since bx^ is arbitrary at any chosen point and F is arbitrary, we have 
the following differential equations for equilibrium: 

(16-7) rS' + pM^ = 0. 

Consequently the virtual work of all the forces (mass as well as surface) 
acting upon any portion of the medium in any virtual displacement is 
given (on using 16-5 and 16-7) by 

(16-8) Total virtual work = J" S J'T^{8x^) dV. 

Since this must vanish for any rigid virtual displacement, i.e., for 

(5a;p).g + (5 Xq),p = 0, 
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the stress tensor must he symmetric: 

( 10 . 9 ) rpuS ^ 

Hence 16.8 can be written 

(16-10) Total virtual work = dV, 

Within the approximations of the usual approximate elasticity theory, 
the total virtual work may be written 

(16- 11) J'J'J'T^ dV 

V 

since formula 15-31 holds for the approximate theory. But note that 
this is not a legitimate result for the finite deformation theory; formula 
16-10 is f^ne legitimate result for that theory. 

Elastic Potential. 

We shall now turn our attention to the elastic potential and its 
relation to the stress tensor. Let p be the density of the volume element 
dV in the strained medium. The element of mass dm is given then by 
dm = p dV. The princi'ple of conservation of mass in a virtual dis- 
placement is expressed by 

5(dm) = 5(p dV) = 0. 

Let T be the temperature of the element of mass dm, <x the entropy 
density (per unit mass) so that the entropy of the mass dm is a dm 
= p<r dV, and u dm the internal energy of the mass dm. Then the 
fundamental energy-conservation law of thermodynamics says that 

(16*12) T'5(<r dm) ~ d(u dm) 

(virtual work of all forces acting on dm). Let 

(16.13) =u - T<r, 

the free energy density or elastic potential. 

From the principle of conservation of mass, we have, on integrating 
over any portion of the strained medium and making use of equations 
16-8, 16-12, and 16-13, 

(16-14) fSJ'{8<t>)p dV = dY - SfS(.^T)p<r dV. 

V V V 

Since Y is arbitrary, this yields 
(16- 15) p 50 = 

We shall now work under the following 

Hypothesis om Elastic Potential 0: 4> is a function of the 
Euclidean metric tensor Qijix) in the strained medium, the Euclidean metric 
tensor ct^c{a) in the unstrained medium, and the temperature T. 
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We shall restrict ourselves to isothermal variations, so that T is a 
constant parameter in From 16 • 15 and the symmetry of the stress 
tensor we see that 50 = 0 for any isothermal rigid virtual displace- 
merd. Now, for any virtual displacement, 5«^c = 0 and dgra = 0. Hence 
from Killing's differential equation 15 • 33 we have 

whenever 

(16*17) -h (5a;^).« = 0. 

Let 

and use the formula 


and 16*17 in 16*16 to obtain 


Z)d> 

^aP 

Lb(«a.p) 


b<f> 


(see formula 15 • 28) 


= 0 . 


Hence <f> must satisfy the following system of partial differential 
equations 

(16*18) 










This is a complete system of three linear first-order partial differential 
equations in the nine variables “a There are nine conditions in 16 * 18 
but three are identities and only three of the remaining six are inde- 
pendent. From the theory of such systems of differential equations ‘ 
we know that the general solution of 16 -18 is a function of six functionally 
independent solviions. It will take us too far afield to give the theory 
of such differential equations ; we are content here with this mere state- 
ment of the result concerning the most general solution 0. 

There are some particularly interesting solutions of equations 16 • 18. 
To consider them it is convenient to define an isotropic medium. 

Definition of Isotropic Medium. A medium whose elastic potential 
is a strain invariant that may depend parametrically on the temperature 
T will he called an isotropic medium. 

Now it can be shown, but in this brief volume we have not the time 
to give the details of proof, that the elastic potential for an isotropic 
medium satisfies the differential equations 16 • 18. It can also be shown 
by a long mathematical argument that any strain invariant is a function 
of the three fundamental strain invariants h, /a, and Zs of Chapter 15. 
The following important result is immediate. A necessary and sufficient 
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condition thcd a medium he isotropic is that its elastic potential 4> = 

I 2 , 1 3 , T) where Ji, J 2 , o/nd Iz ore the fundamental strain invariants. 

In the usu al approximate theory of elasticity, the elastic potential 0 
for crystalline media (another name for non-isotropic media) is taken 
as a quadratic function of the strain tensor components. It is tacitly 
assumed in the usual approximate theory that a special privileged 
reference frame, determined by the axes of the crystal, has been 
chosen. The coefficients of the quadratic form are accordingly not 
scalars but components of a tensor of rank four which depends on the 
orientation of the crystalline axes. 


Stress-Strain Relations for an Isotropic Medium. 

Consider the elastic potential 0 for an isotropic medium as a function 
of the strain tensor €ra. Since is symmetric, we have = !(€« + €.r). 
In we shall write ^(era + u-f) wherever era occurs, and thus we see that 


(16-19) 


' ' ' " = — ‘ — j 

^€ra ^^ar 


with the understanding that in - — , say, all the other e^s (including 

^€7*5 

for that s, r) are held constant, so that in this differentiation no atten- 
tion is paid to the symmetry relations €ar — €ra- 

We saw in Chapter 15 (see 15-29) that under a virtual displacement 
the variation of hp^ and hence of the strain tensor epg was given by 

(16-20) dcpg — —2^hpq = 


since hpq = Qpg — 2€pq. But 5gpg = 0; hence 

50 = ~ Oip i5Xr),cc + hi (Sa:,.)^], 

or 


(16-21) d<i> 


b0 




(a and /3 are summation indices t) 


on using conditions 16 • 19. Now, for an isothermal virtual displace- 
ment, formula 16-15 reduces to p 50 = T‘“^(5a;„),^, and so for an isotropic 
medium 


(16-22) 


2)0 






From the arbitrariness of the virtual displacement and the fact that 


t Throughout the remaining part of this chapter, a mere repetition of an index 
in a term will denote summation over that index. 



94 


STRESS TENSOR, ELASTIC POTENTIAL 


= 6 ? - 
(16*23) 


2eJ we obtain the stress-strain relations for an isotropic medium 



We put this stress-strain relation in another form. Now by 
definition ^ and hence 



II 

/O I/O 

Obviously 

b<j> b4> bej; 

be^- beij 

and hence 


(16*24) 

b0 w 

ben ' 


Define the mixed stress tensor T% by 
(16*25) 

(= from the symmetry of the stress tensor). Then with the aid 

of 16 • 24 in 16 • 23 one can show readily that the following stress-strain 
relations hold for an isotropic medium 


(16*26) 






From the principle of conservation of mass p dV — po dFo, and from 
the fundamental result 15 * 5, we see that 

(16-27) p = poVl - 2h + 4/2 - 8h 


in terms of the strain invariants Ii, 72, and /.j for media, whether iso- 
tropic or not. Since 7i, I 2 , and are respectively first degree, second 
degree, and third degree in the strain tensor comi^onents we see 
that to a first approximation p = pu; i.e., volumes are also preserved to 
a first approximation. Hence to the same degree of approximation, the 
stress-strain relations 16-26 for an isotropic medium reduce to Hookers 
law of the usual approximate theory 



(16*28) 

where = p0. 
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TENSOR CALCULUS IN RIEMANNIAN SPACES AND THE 
FUNDAMENTALS OF CLASSICAL MECHANICS 


Multidimensional Euclidean Spaces. 

In the last two chapters of this book we shall attempt to give some 
indications of a more general tensor calculus and some of its applica- 
tions. Although our discussion will of necessity be brief, this fact will 
not keep us from going to the heart of our subject. Our study of 
Euclidean tensor analysis can be used advantageously to accomplish 
this. 

First of all the subject matter of Chapters 9, 10, and 11 can obviously 
be extended to rt-dimensional Euclidean spaces, where n is any positive 
> integer. There will be n variables wherever there were three before, 
and indices will have the range 1, 2, • • • to n with the consequent 
summations going from 1 to n. For example, the squared element of 
arc in rectangular coordinates 2/\ • • • , 2/” is 

(17*1) 

i— 1 


■while in general coordinates x'-, x^, - ■ a:" 


(17-2) 

where 

(17-3) 


ds* = gaftix^, x^, • • • , a:’*) 


gaf,ix\ x\ s») 


h dx" dx^' 


the re-dimensional Euclidean metric tensor (see 9- 15). The 7r-dimen- 
sional Euclidean Christoffel symbols are 


(17-4) 


ri^(a:S x‘, 




a<r 



where the ffa/s Q-re defined in 17*3 while 


(17-5) 
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in terms of the w-rowed determinant 


( 17 - 6 ) 


gii, gizj 
g2i, goz, 


•> ffln 

‘r g^n 


gnlj gn2f * * *> gnn 


Geometry.^ 

An n-dimensional Riemannian space is an n-dimensional manifold 
with coordinates such that length of curves is determined by means of 
a symmetric covariant tensor field of rank two ga$i^^y 
such a fashion that the squared element of arc 

( 17 - 7 ) 

is positive definite, i.e,, g^^ dx^ dx^ ^ 0 and is equal to zero if and only 

The length of a curve = p{t) given in terms 
of a parameter t is then by definition 


if all the dx"* are zero. 



sc* (*i) 


Fia. 17-1. 


( 17 - 9 ) 


g 


It can be proved by rather long alge- 
braic manipulations that from the positive 
definiteness of g^p dx^ dx^ follows the posi- 
tive value of the determinant of the g^fi, i.e., 

g^ii, S'12, * * * , gin 


> 0 . 


gnl] gnZ} * * * j gnn 


The theory of a Riemannian space is a Riemannian geometry. 

Exactly as in an n-dimensional Euclidean space, we can derive the 
contra variant tensor field of rank two o?, • - • , x'^') and thus have 

at our disposal the Christoffel symbols of our Riemannian geometry 


( 17 * 10 ) 



(a5^, x^, - • - , X") 


1 ia (^9^ . 

2 'bafi bx'^ / 


Notice that the Riemannian Christoffel symbols depend on the funda- 
mental Riemannian metric tensor g^^ and its first partial derivatives. 

Unlike the Euclidean Christoffel symbols, it is impossible in general 
to find a coordinate system in which all the Riemannian Christoffel 
symbols are zero everywhere in the Riemannian space. This is due to 
the fact that it is in general impossible to find a coordinate system 
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in which the fundamental metric tensor has constant components 
throughout space. It is to be recalled that in a Euclidean space there 
do exist just such coordinate systems, i.e., cartesian coordinate systems 
and rectangular coordinate systems in particular. 

We can, however, prove that there exists a coordinate system with 
any point of the space as the origin, i.e., the (0, 0, - • • , 0) point, such 
that all the Christoffel symbols vanish at the origin when they are 
evaluated in this coordinate system. Such a coordinate system is called 
a geodesic coordinate system. We shall prove that the coordinates 
defined implicitly by the transformation of coordinates 

(17-11) + • • - , 


are geodesic coordinates. 

A direct calculation from 17-11 yields the needed formulas 



where the 0 means evaluation at the origin of the y^ coordinates. Now, 
under a transformation of coordinates, the Christoff el symbols of a 
Riemannian space transform f by the rule 10-29 for Euclidean Chris- 
toffel symbols. Let y^, • • * , y") be the (Riemannian) Chris- 

toff el symbols in the y^ coordinates. Then 


(17-13) *Vip{y\ 2/^ * • *, T) = ri„(a;S x\ - - - , a;«) 




b^“ ba;^ 


Now we see from the transformation of coordinates 17-11 that a;* = g* 
when 2/^ = 0. In other words, the origin of the y^ coordinates has co- 
ordinates X'- - in the x'- coordinates. 

If we evaluate both sides of 17 - 13 at the origin of the 2/* coordinates 
and if we use formulas 17-12 in the calculations, we find 

(17-14) 

In other words, the 2 /^'s are geodesic coordinates. 


t With the difference that the number of variables now is n and the indices 
have the range 1 to n. The proof is practically a repetition of that given in note 4 
to Chapter 10. 
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Curved Surfaces as Examples of Riemannian Spaces. 

Obviously any Euclidean space is a very special Riemannian space. 
A simple example of a Riemannian space which is not Euclidean is 
furnished by a curved surface in ordinary three-dimensional Euclidean 
space. This can be seen as follows. Let ?/* be rectangular coordinates 
in the three-dimensional Euclidean space, and let the equations of a 
curved surface be 


(17*15) 

in terms of two parameters and x^. Then the squared element of arc 
for points on the surface 17 • 15 is 

(17 • 16) (d2/0^ = dx“ daf 

i=l 


(a and /S have the range 1 to 2 and corresponding summations go from 
1 to 2 ), where 


(17 • 17) 




^ b/^(3;S x^) bP(x^, x^) 


So a surface in three-dimensional Euclidean space is a two-dimensional 
Riemannian space. 


Exercise 


The surface of a sphere is a two-dimensional Riemannian space. 
Find its fundamental metric tensor and its Christoffcl symbols. The 

surface of a sphere of fixed radius r is 
given by 

— r sin x^ cos x^ 
ip- = r sin sin x^ 
ip = r cos 

Therefore the fundamental metric tensor 
is given by gn = = ^21 = 0 , 

9^22 = r2(sin x^p. Hence 

The Christoffel symbols are then 

TI 2 = —sin x^ cos x^, rfa = Pai = cot x^, 

and all the other Christoffel symbols of the surface of the sphere are 
zero. 
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The Riemann-ChristofEel Ctirvature Tensor. 

It was seen in Chapter 1 1 that covariant differentiation is a commu- 
tative operation in three-dimensional Euclidean space, and by exactly 
the same type of reasoning this is also true in an n-dimensional Eu- 
clidean space. To establish this result explicit use was made of car- 
tesian coordinates. Since such coordinates are in general not available 
in a Riemannian space, we cannot use that type of proof. In fact, co- 
variant differentiation in a Hiemannian space is not in general commu- 
tative. We shall find a formula (see 17 • 19 below) that makes clear the 
non-commutativity of covariant differentiation in Riemannian spaces. 

In obtaining Laplace’s equation in curvilinear coordinates for con- 
travariant vector fields in a Euclidean space, we had to calculate the 
second covariant derivative of a contra variant vector field. (See the 
bracket term in 12*23.) The calculation for Riemannian spaces is 
practically the same, so that we shall write down the second covariant 
derivative of x^, • • •, a:") based on the (Riemannian) Christoff el 

symbols r*^(a;^, rc®, • • *, x'^') without giving any more details (again 
see bracket term in 12*23). The result is 


(17*18) 


da:“ 




I T-ll 


51 


\ 


+ r 


X* -pT 

TjS**- act 


r % pr 
crr^ oe/3 




From the commutativity of the partial derivatives and the symmetry 
of the Christoffel symbols we find 

(17-19) 


where 

(17*20) 



br* 

o-Qg 

Z>x^ 


dr* 


<r/3 


da:" 


. p’- 






To justify the notation and prove that they are the com- 

ponents of a tensor field of rank four, contravariant of rank one and 
covariant of rank three, we first note that the left sides of 17 • 19 are 
the components of a tensor field of rank three, contravariant of rank 
one and covariant of rank two. Hence is a tensor field of the 

same type for all contravariant vector fields i.e., 

_ . Z)x'' d:r* 


But, writing $^(a;) in terms of ^(x), we evidently have 




do:'* dx" dx" d^* 
dx"^ dx“ dx^ 
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But |‘'(») are arbitrary, and hence, equating corresponding coefficients, 
we obtain the tensor law of transformation for This tensor 

field is the famous Riemann-Christoffel curvature tensor; it is not a 
zero tensor in a general Riemannian space. Hence ^ ® 

Riemannian space with non-vanishing Riemann-Christoffel curvature 
tensor. But obviously the Riemann-Christoffel curvature tensor is 
zero in Euclidean spaces. Hence Euclidean spaces; 

this checks a result found earlier, in 11.13. 


Geodesics. 


A straight line is the shortest distance between two points in Eu- 
clidean spaces. There are curves in Riemannian spaces that play a 
role analogous to the straight lines of Euclidean spaces. Such curves 
are called geodesics. In fact, if a Riemannian space is a Euclidean 
space, then its geodesics are straight lines. To find the differential 
equations satisfied by the geodesics of a Riemannian space, we have 
to get the Euler-Lagrange differential equations for the calculus of 
variations problem 


(17-21) 


j dx* dx^ 

A \ It 


dt 


minimum. 


It can be shown that the Euler-Lagrange equations for this calculus 
of variations problem are 


(17-22) 


cPxKs) 




dx^^dx^ 
ds ds 


where 5 is the arc length and r«^(x) are the Christoffel symbols of the 
Riemannian space. In other words, if the coordinates of points on a 
geodesic are considered as functions x^{s) of the arc length parameter 
s, then the n functions x*(s) satisfy the system 17 • 22 of n differential 
equations of the second order. 

If the Riemannian space is Euclidean and we choose rectangular 
cartesian coordinates y^, equations 17 • 22 reduce to 


and hence 


d^*(g) 

ds^ 


y* = oi^s + /3*' (a* and jS® are constants), 

the parametric equations of straight lines in terms of arc length s. 
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Equations of Motion of a Dynamical System with n Degrees of 
Freedom. 

In classical meclianics, it is postulated that the motion of a con- 
servative dynamical system of u degrees of freedom with no moving 
constraints is governed by LagTange^ s eguations of motion^ 


(17-23) 



If the kinetic energy T, in terms of the generalized coordinates g\ 
• • •, g” and the generalized velocities g‘(0 = ^ j is 

T = ^gii(gS • • •, (Qa = 


and if the potential energy is V (gS g*, • - • , g”) , then the kinetic potential 
or Lagrangean L is given by L = T — V. Now the kinetic energy is 
positive definite in the velocities g*j i.e., T ^ 0 and T = 0 if and only 
if gi _ Q If ean be proved by algebraic reasoning that the determinant 
g of the gij is positive so that we can form in terms of the Qa ex- 
actly as in Riemannian geometry. By direct calculation we find 



Hence Lagrange^s equations of motion 17*23 can be written in the form 

(17-24) ffi// &?• 

Multiplying corresponding sides of 17-24 by g“* and sunaming on i, 
we obtain the following form for Lagrcingo^ s eguation of motion^ 


(17-25) g“ + r^/c(gS * • • , q^Wq^ = 

where r“A;(g^ - - • , g”) are the Christofel symbols based on the gij (gS 
* * * j 9”) kinetic C'nergy of the dynamical system. 
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Exerdse 


A symmetrical gyroscope with a poiat 0 fixed on the axis is acted upon by grav- 
ity. Let I, I, andJ be the principal moments of inertia. Then the kinetic energy 
is given by 


T.-JP- 
2 


+^/sin^j' 



1 /df 


df^ 

1 } 


and the potential energy by 

V = Mgh cos jK The coordinates q\ q\ and q^ are the Eulerian angles, M is the 
mass of the gyroscope, and h is the distance of the center of gravity from 0. Find 
the Lagrangean equations of motion of the symmetrical gyroscope. Compute also 
the element of arc length of the three-dimensional Riemannian space associated 
with the symmetrical gyroscope. 



CHAPTER 18 

APPLICATIONS OF THE TENSOR CALCULUS TO 
BOUNDARY-LAYER THEORY 


Incompressible and Compressible Fluids. 

The constancy of volume of all parts of a fluid in motion sometimes 
plays an important role in tlie theory of fluid flows. A fluid in motion 
with this property is called an incompressible jiuid, whereas a fluid in 
motion without this property is called a compressible fluid. If w* are the 
contravariant components of velocity of the fluid in motion in general 
coordinates x*, then 


(18.1) 



are the differential equations whose integration gives the paths of the 
fluid particles in the coordinates x\ 

It can be proved by a direct calculation that a necessary and suflflcient 
condition that the volume 

(18 . 2) S dx^ dx^ 

of arbitrary portions of the moving fluid be preserved ^ is that the 
divergence of the velocity field be zero, i.e., 


(18-3) < = 0. 

In 18*2, g is the determinant of the Euclidean metric tensor 
(ds^ = gij dx^ dx^'), and the comma in stands for covariant differen- 
tiation based on the Euclidean Christoffel symbols r}/^. In other words, 
CL necessary aiid su jjicient condition for an incoynpressible fluid is thcd the 
velocity vector field u^ satisfy the partial differential equation 18-3. A 
glance at formula 13 • 6 shows that the condition of incompressibility is 
equivalent to 

(18-4) = 0. 


If we recall the Navier-Stokes equations 13*3 for the motion of a 
viscous fluid, incompressible or compressible, we know that the equation 
of continuity 


(18-5) 
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in general coordinates a: » merely states the constancy^ of the mass m 

(18-6) m = S S S dx^ dx^ 

of any portion of the moving fluid. An evident consequence of the con- 
ditions 18-3 and 18*5 is that the density p(jx^, a^, t) (an absolute 
scalar) satisfies the condition 


(18*7) 

which states that 
(18*8) 


dp dp „ _ 


dp 

dt 


= 0 


along any chosen path of fluid particles. This means that 18-7 can he 
taken as the defining condition for an incompressible fluid in view of the 
continuity equation 18.5. The Navier-Stokes equations for an incom- 
pressible viscous fluid reduce then to the following system of four 
differential equations in general coordinates a;*: 


(18*9) 


{ 


du" 

d^ 

a 


= 0. 


u^u\ 


1 XT- 

7 5 ? + 


For an incompressible fluid, the density p(a^, x^, a;®, t) is given subject 
to condition 18 • 7. Then the four differential equations 18 • 9 will have 
as unknowns the three velocity components w®, w® and the pressure 
p(pc^, x^, a®, t) of the fluid. 

The situation is different for compressible fluids. The Navier-Stokes 
equations 13.3 are four in number with five unknown functions u^, w®, 
w®, p, and p. To make the problem determinate a fifth condition must 
be imposed. This is usually furnished by the “equation of state, ’’ 
which in the isothermal case is of the form 


(18*10) 


V = /(p)* 


Boundary-Layer Equations for the Steady Motion of a Homogeneous 
Incompressible Fluid, f 

We shall now restrict ourselves to the steady motion of a fluid with- 
out any external forces, so that A* == 0 and all the quantities u\ p, p are 
independent of the time t. If in addition we assume that the fluid is 
homogeneous, i.e., p is a constant, and incompressible, the four unknowns 

t The remaining part of this chapter is an exposition of some unpublished re- 
searches of Dr. C. C. Lin. These results were presented by Dr. Lin in my seminar 
on applied mathematics. 
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w*, p must, by a reference to 18 • 9, satisfy the four differential equations 


(18*11) 


{ 


- g 

= 0 , 





where x is the pressure p divided by the constant density p of the 
fluid, and the constant v is the kinematical viscosity. Since the co- 
variant derivative of the Euclidean metric tensor gij is zero, it follows 
from 18*11 that the covariant vector components Ui = QiaU* of the 
velocity field and the function x will satisfy the system of differential 
equations 

( „ &x 

(18-12) ^ 

( w, = 0. 

For the treatment of ‘^boundary-layer’’ problems connected with 
an arbitrary surface, it is convenient 
to take a system of space coord/mates 
in which x^, x^ are surface coordinates 
and a? is a coordinate measured along 
the normals to the surface. Thus 
= 0 will be the equation of the given 
surface. If we allow Latin indices to 
run over the range (1, 2, 3) and Greek 
indices over the range (1, 2), we have 
the following fundamental metrics: 

(18*13) cc®) dx^ dx^ in 3-space, 

and 

(18*14) = Sf«^(xS ^2, C) dx^ dxP 



over the surface 

a;® = C, a constant. 

From the manner in which the coordinate x? was chosen, it follows that 
in 3-space t 

(18*15) ds® = gij dx^ dx^ = + {dx^Y, 


t In Riemannian geometry, this is sometimes called the **geodesic form” of the 
line element. Such forms of the line element were used recently by Dr. W. Z. Chien 
in connection with his researches on the intrinsic theory of plates and shells (see 
references). Dr. Chien presented some of his work in my seminar on applied mathe- 
matics and made some exceedingly helpful calculations in connection with interest- 
ing geometric ideas arising in his and Dr. C. C. Lin’s work. 
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t 

In other words, the Euclidean metric tensor x^, x?) is such that 

( 18 * 16 ) gzz = 1 , gz^ = ^«3 = 0 . 

We shall henceforth consider transformations of surface coordinates x^ 
alone so that the coordinate x^ may be regarded as a scalar parameter 
under a transformation of surface coordinates. To emphasize this 
fact we shall use the notation x° = Xq — x^. Thus, in the new notation, 
18 • 16 can be written 


(18 • 17) ^00 = 1, ^o« = gao = 0. 

We saw in the previous chapter that a surface can be considered as a 
two-dimensional Riemannian space. There is thus at our disposal the 
Riemannian tensor calculus of the previous chapter for immediate use 
in connection with the surface re® = constant. We shall use a semicolon 
to denote surface covariant differentiation in contradistinction to the 
comma for covariant differentiation in the enveloping three-dimensional 
Euclidean space. 

To express all covariant differentiations with respect to space co- 
ordinates in terms of covariant differentiations with respect to surface 
coordinates x^, x? and partial differentiations with respect to rr®, con- 
sider first the Euclidean Christoffel symbols Py* in the coordinates 
x^, x^, X?. If i, j, k are all in the range 1, 2 no reduction is possible 
unless a special surface coordinate system is chosen; if one of the 
three indices is zero, we have 


(18-18) 


•pO 

J- 


2 dxo ’ 


Tia 

— A/30 


— ^ 


^0y 


These are evidently tensor fields with respect to transformation of 
surface coordinates, and they shall be denoted by and r| respec- 
tively. The other Christoffel symbols P}*, in which two or all of the 
ij jf k are zero, vanish identically. 

With the help of these relations, it can be easily verified that 


(18-19) 

(18-20) 


Won — 




2>Wo 


ba;" 




~ “ ^«;/8 + TaisWo, 




(18-21) 














4-^. 


ca 
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where 


! 


«>o = + 2r“V:fl - rirSuo + r^u„ 


(18-22) = ^TlV%u^ +riu^ + g^'^Ucc-,p;y “ - Ti,^Uo 

- (rjr? + r^r2)«s, 


do not involve differentiation of Ui with respect to a^. 

Let us now consider the analytical nature of the system 18-12 of 
four partial differential equations in the four unknowns t, Uq, By 
using 18-19, 18-20, and 18-21, we can put this system in the normal 

form with respect to aP by solving for from the equations 

of motion and for from the equation of continuity. Thus 




(18-23) 


bTT 

bx^ 


^Uq _ 

bx^ “ bx^ ^ 


a®" 




bTT 

b^ 


ft 

r|r-? + 4‘, 




buo 


= + r^'Wo, 


where and in the first equation may be expressed in terms of 


bx^ 

bu. 


bx^^ 


Uq and r— ^ by using the last equation. Thus, the highest derivatives 
Ox^ 

of all the variables with respect to have the coefficient unity in 

these equations. Hence, if ir, Wo, are given as functions of 

and x^ on the surface x^ = 0, the solution of the problem is imiquely 
determined. 

This normal form 18-23 of the system of differential equations, how- 
ever is not analytic in the small parameter x', the important case in 
aeronautics, in the neighborhood of x' = 0, and is consequently incon- 
venient for the application of the method of successive approximations. 
We therefore make the transformation of variables 


(18-24) 






■\/ V 


w — 


Uq 
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to bring it into the desired form. We then have 

^ ^ ft / a/ 


5pr-r»5p + 4>o^ 


(18-25) 


br \ot Ox“/ W < 

_ ^ + „v-.a + - 2r£ — 

. r«— 

^ = -wf^ + j^r^. 


Let us note that and <l>« are linear in the small parameter through 
the term in uo (of. 18 • 22), and also depend on through the geometrical 
quantities, which are functions of x° — Indeed, it can be shown ^ 
that the surface metric tensor is a quadratic function of x^, while all 
other geometrical quantities may he expanded as power series of xP con- 
vergent for I < Rmj Rm being the minimum magnitude of the 

principal radii of curvature over the surface under consideration. 
Hence, the right-hand sides of the equations 18 * 25 are Taylor series in 
v^, and the solution of 18 • 25 may be carried out by expanding each of 
the dependent variables as a power series of v^, convergent for all finite 
values of for which | | < Rm. 

If we try to solve 18*23 by the same type of expansion, either the 
series are asymptotic, or they may terminate; but in general wo cannot 
find a solution satisfying all the required boundary conditions. In 
fact, the initial approximation is easily verified to satisfy the non- 
viscous equation (r = 0). The boundary conditions at infinity and the 
condition Wo = 0 at = 0 are then sufficient to determine tins approxi- 
mation completely. Indeed, the boundary conditions at infinity are 
usually such that the resultant solution is potential. Then the initial 
approximation is an exact solution of the complete equations 18*23. 
However, the boundary conditions Wo = 0 at xP = 0 cannot be satisfied 
in general. The effect of viscosity can never l.)e l)rought into evidence. 
This shows that the more elaborate treatment described above is absolutely 
necessary. The non-viscous solution (usually potential), however, 
serves as a guide for making the exact solutions satisfy the boundary 
conditions at infinity. This point will be discussed in more detail 
below. 

Let us now proceed with the solution of 18.25 by writing 


i TT — + VTT^^'^ + * * • + * * *, 

< Wa = -h + * * * -1- - ■ - , 

i^ w == PW^^^ + * • • -j- • • - . 


( 18 - 26 ) 
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Corresponding developments for the geometrical quantities 
must also be used. The initial approximation gives 


(18-27) 


+ V) 




0 

Zyw 


Stt 

= 0 , 


+ 


br' 


where the superscripts of the initial approximation are dropped. In 
these equations, T is a scalar, and the metric tensor is being 

Qa^^y ^y cvaluatcd at rr® = 0. The conditions over the surface 
^ = 0 are Wq, = 0 and ly = 0. The condition at infinite f is set according 
to the following considerations. For a large but finite value of 
the value of ^ is still small. Hence, the solution may be expected to 
pass into the non- viscous solution close to the surface if f is large. 
Thus, for the initial approximation, we may lay down the conditions 


(18 • 28) Ua, = TT = -K for f CO , 

where and tt are functions of and x^, being the values of and tt 
of the non-viscous solution at x^ = 0. The initial approximation is 
then completely determined. 

If Uc and TT differ from and x by quantities of the order of v for 
t = h, then an approximate solution of 18- 12 is usually taken to be 
given (a) by the non-viscous solution for f > ^, and (h) by the solu- 
tion of 18-27 for ^ <. h. The quantity li^/ v is known as the ^‘thickness 
of the boundary layer” and is arbitrary to a certain extent. For ex- 
ample, we may define h to be given by (say) three times h of the 
equation 

J f CO 

- uj dt = ujiy 

0 


which is in general different according to whether a = 1 or 2. This 
initial approximation is usually known as the houndary-layer theory of 
Prandtl. Incidentally, we note that tt is a function of and x^ alone, 
by the second equation of 18-27. Hence, by 18-28, tt = W ( x '‘-, x ^), 
which is known from the non-viscous solution. The first and third 
equations of 18-27 then serve as three equations for the velocities 
Ua and w. 

The higher approximations in 18 - 26 satisfy certain differential equa- 
tions obtained together with the derivation of 18-27. The boundary 
conditions at f = 0 are u^ — 0, for any approximation. The boundary 
conditions for the nth approximation at infinity will be specified by 
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using the nth approximation of the asymptotic solution, which will in 
turn be determined from certain boundary conditions related to the 
{n — l)st approximation of the convergent solution. Since we are 
never concerned with higher approximations in practice, we shall not 
go into further details. 
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Chapter 1 

1. In the modem quantum mechanics of theoretical physics, matrices with an 
infinite number of elements as well as with a finite number of elements are used 
very widely. The elements of these matrices are often complex numbers. The 
reader is referred to the bibliographical entries under Born, Jordan, and Dirac. 

2. For applications of matrices to the social sciences the reader is referred to 
the references given in a recent paper by Hotelling. 

3. We shall deal for the most part with matrices whose elements are real or 
complex numbers. It is possible, however, to deal with matrices whose elements are 
themselves mairices. W^e shall have occasion to use a few such matidces in connection 
with our discussion of aircraft flutter in Chapter 7. 


Chapter 2 

1. For the properties of determinants, linear equations, and related questions 
on the algebra of matrices, see Bocher's IntrodxaMon to Higher Algebra. 

2. Cramer’s rule for the solution of linear algebraic equations is given in most 
books on algebra. In our notations it can be stated in the following manner. If the 
determinant u => | aj \ of the n eguations 

ay = ¥ 

in the n unknowns x\ s?, is not zero, then the eguations have a unique solution 

given by 

A' 

X* = — » 

a 

where A* is the n-rowed determinant obtained from a by replacing the elements oj, af, 
•■■,a^ of the ilh column by the corresponding elements 6S 6*, • * *, 6". 

3. The rule for the multiplication of two determinants takes the following form 
in our notations. // a = [ aj [ and 6 « [ 6j 1 are two ur^owed determinants, then 

the numerical product c ab is ttself an n~rowed determinant with elemsnts Cj given 
by the formula 

c) - aibf. 

3 ct 3 

4. It can be shown that Sr, the trace of the matrix A^, is also equal to the sum 
of the rth powers of the n characteristic roots of the matrix A. 

5. The recurrence formula 2*6 for the coefficients ai, ■ • • , an of the characteristic 
function of a matrix can be derived from Newton’s formulas; see Bocher’s Introduo- 
turn to Higher Algebra, pp. 243-244, for a derivation of Newton’s formulas. 

Ill 
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Chapter 3 


1. To prove that the matric exponential is convergent for all square matrices 
A, let A = II aj || be an n-rowed square matrix, and let F(A) be the greatest of the 
numerical values of the n* numbers aj — the greatest of the moduli of the ai if the 

aj are complex numbers. Then each element in the matrix A * will not exceed n*~^V* 

in numerical value. Hence each of the n* infinite series in wiU be dominated by 
series 


, „ nV^ n®F3 

1 4. P' -I _f. 4_ 

X -I- r -r 2 ! -t- 3 J 


n 


Hence all the n* numerical series in converge. This means that is convergent 
for aU square matrices A . 

In the terminology of modem functional analysis and topological spaces, the 
F(A) is called the norm of the matrix A, and the class of n-rowed matrices with 
the operations of addition and multiplication of matrices, multiplication by num- 
bers, and convergence of matrices defined by means of the norm F( A ) — in other 
words, F(A) plays an analogous role to the absolute value or modulus of a number 
in the convergence of numbers — is called a normed linear ring. Other equivalent 
definitions of the norm of a matrix are possible. For example, F(A) can be 
taken as 



if the elements o*. of the matrix A are real numbers. Whatever suitable definition 
of a norm is adopted, the norm F(A) of a matrix will have the following properties: 

(1) F(A) ^ 0 and = 0 if and only if A is the zero matrix. 

(2) F(A + B) ^ F(A) + V(B) (triangular inequality). 

(3) V(AB) V(A)V(B). 

From property 3 it foUows that F(A") (F(A))", a result that makes obvious the 

usefulness of the notion of a norm for matrices in the treatment of convergence 
properties of matrices. 

The class of matrices discussed above is only one example of a normed linear 
ring. The first general theory of normed linear rings was initiated in 1932 by 
Michal and Martin in a paper entitled “Some Expansions in Vector Space,” Journal 
de TnalhArnaiiqueB pures et appliquSes. 


2. The special case of the expansion 3 • 7 when F(A ) is a matric polynomial is 
known as Sylvester’s theorem. 

If the characteristic equation of a matrix Ahas multiple roots, then the expansion 
3-7 is not valid. However, a more general result can be proved. For the case of 
matrix polynomials F(A) see the Duncan, Frazer, and Collar book. The more 
general cases of matric power series expansions arc treated briefly in a paper by 
L. Fantappie with the aid of the theory of functionals, since the elements FJ(A) 
in the F(A) bxg functionals of the numerical function F(X). See Voltcrra's book on 
functionals (Blackie, 1930). 


3. Another equivalent form for the n matrices Oi, (? 2 , • 

M(Xi) 


•, (S'n is 




/d/(X)\ 

\ dX /\=\i 
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where the X» are the characteristic roots of the matrix A, /(X) is the characteristic 
determinant of A, and M(\) = || ilfy(X) |1 is a matrix whose element is the 

cofactor of —• in the characteristic determinant /(X). Notice carefully the 
position of the indices i and j. 


Chapter 4 

1. A good approximation to the solution X(t) = of = AX(t) 

at 


(t — 

can be obtained by taking 1 + (t — ta)A + — — — + 


(< — <o)" , . 

•i : — A”- m the 


2! ■■ ■ ■ n! 

place of the infinite expansion for Pqj, many practical purposes n = 2 would 

be large enough to give a good approximate solution. The approximate solution 
can then be written 

Xit) ^Xo + it- U)AXo + - ~l —A^Xc, 
where X^ is the column matrix for i = io initially given. 


2. If the solution of the matric differential equation 

(1) = AXif) (XOo) = Xo) 

has been found, then the solution of 

(2) ~U+ bI)X(.t) (X(to) = Xo) 

at 

can be written down immediately. In fact, from the second property of the matric 
exponential given in Chapter 2, we see that where e* is the numerical 

exponential. Hence by formula 4 • 3 we see that the solution of equation 2 is ob- 
tained by a mere multiplication by e** of the solution of equation 1. 


Chapter 6 

1. The reader is referred to Whittaker’s Analytical Dynamics for a treatment 
of Lagrange’s differential equations of motion of particle dynamics. For some 
engineering applications, the reader is referred to Mathematical Methods in Engi- 
neering by Kdrmdn and Biot. 

2. Consult the references in the above note. 


Chapter 6 

1. If only the fundamental frequency is wanted but not the corresponding 
amplitudes, then the application of Rayleigh’s principle may be preferable. See 
Elementary Matrices by Frazer, Duncan, and Collar, pp. 310, 299-301. 

2. A special case of Sylvester’s theorem is what is actually used; see expansion 

3*9. 
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Chapter 9 

1. Although little use has been made of the tensor calculus in plastic deforma- 
tions, one would suspect that a thoroughgoing application of the tensor calculus 
to the fundamentals of plastic deformation theory would prove fruitful. 

2. For some elementary applications of the tensor calculus to dynamic meteor- 
ology, the reader is referred to Ertel’s monograph. 

3. We shall deal briefly with Riemannian spaces (certain curved spaces) and 
their applications to classical dynamical systems with a finite number of degrees 
of freedom (see Chapter 17), and to fluid mechanics (see Chapter 18). 

4. A discussion of the fundamentals of coordinates, coordinate systems, and 
the transformation of coordinates in the various spaces, including Euclidean spaces, 
is out of the question here. The readers who are interested in modem diflerential 
geometry and topology will find ample references in the bibliography under the 
entries for Veblen, Whitehead, Thomas, and Michal. 


Chapter 10 

1. Some writers, especially those dealing with physical applications, like to think 
of the contra variant and co variant components of one object called a vector. For 
example, if are the contra variant components of a velocity vector field, then 
and can be considered the contravariant vector and co variant vector “repre- 
sentations” respectively of the same physical object called “velocity vector field.” 
This point of view, however, is untenable in spaces without a metric gij. 

2. The importance of the Euclidean Christoffel symbols for Euclidean spaces 
is, even now, not very well known. 


3. Since 

a®, ««) 




(l/S are rectangular coordinates). 


we obtain, from the rule for the multiplication of two determinants, the result that 

g = I gap I = where J = 

determinant, of the transformation of coordinates to rectangular coordinates y* 
from general coordinates x*. Hence J 0 since we deal with transformations of 
coordinates that have inverses. This means that the determinant p 0 for all our 
“admissible” transformations of coordinates. 






is the Jacobian determinant, or the functional 


4. The following steps establish the law of transformation 10*29 of the Eu- 
clidean Christoffel symbols; exactly the same method establishes the corresponding 
law of transformation for the Riemannian Christoffel symbols to be discussed in 
Chapter 17. 

Since Qixv are the components of the Euclidean metric tensor we have under a 
transformation of coordinates from coordinates x* to coordinates £* 

(«) g<rp(.3i\ 53 ) = gM.Cx\ x^ x^)— 
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Differentiating corresponding sides of (a) we obtain (considering the as inde- 
pendlent variables) 


(b) 

^9 <tP 

b^pv bx** bx" bx^ r b^x^ bx"”] 

bx^* 

b*x" 

b2“ 

“ b^ bi® bi^ bx« [_^^‘’b2® bx« b#_ 

‘'^b2® bx^ b2“ 

Interchange a 

and <r in (6) and, noting that 





b^x^ b^x^ 





bx® bx“ bx“ bx® 



obtain 





(c) 

^QaP 

bgp^y bx^ baf bx'' F b^x'* ba:"”j 

r bx'* 

I 1 

b^x" 1 

b2® 

bx'' bx“ bx^ b2® 5x® bx® b2^J 


bx^ b2®_j 


Interchange /3 and o- in (c) and, noting that 




obtain 

(d) 


^Qjxv baj^ 


bx^ bx^ b2“ bx"” bx^ 


zzff + fl'Mi 


b^x^ bx^ 


bx^ bx“ bx*" 


1^ + 


[ 5x^ b^af T 
^^bSf^bsP b2®’J* 


Add corresponding sides of (6) and (d) and subtract corresponding sides of (c) after 
interchanging fi and p in the first terms of (6) and after interchanging v and p in the 
first terms of (d). Then take -J- of both sides, obtaining 

1 / b^<rg bgcta _ b^\ _ 1 / ^ ^ j_ 1 „ ^ 


(e) 


2Vb2« bx^ bx*"/ 


2 Vbx'" ' bx" bx'* y b2“ bSP b£^ 2 b# b2“ 

1 b^x^ bx^ 
■^2^*^b2<3 52“b2'^’ 


the terms enclosed in brackets in (6), (c) and (d) canceling out in the additions and 
subtractions. On interchanging ju and v in the last terms of (c) and on recalling that 
Ovti = we get 


(/) 
Now 
iff) 


1 / b£^ bgacr _ b£^\ ^ 1 / b^ bflr^p _ b£^\ ^ ^ ^ ^ 


2Vbx“ b2^ bx*^/ 2Vbx'^ 


bx" 


bx^ J bx“ b2^ bx®" ^^’’bx® bx^ b2“ 






,lo- 




,X<«} 


bx* bx® 
bx^ ba^ 


Multiplying corresponding sides of (/) and {g), summing on <r, and using the identi- 
ties 


we readily obtain 


gial 

! Vbx“ bx^ 


bx® / 


bx'* bx® 

=: 5^ 3 

b2®^ 

tiff ^ 



Ibx'* 


^*'ffM» 




MP 


b^pv\ ba:'* bx" b2* b^x^ b2* 


0 


da;'' / d:^ dx^ dx^ dx^ d£“ dx^ 


from which the desired transformation law 10*29 follows immediately on recalling 
the definition of and 


4 
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Chapter 11 

1. Formiala 11-12 for the co variant derivative of a tensor field can be estab- 
lished very quickly with the aid of normal coordinate methods of modern differential 
geometry. See the two Michal and Thomas 1927 papers. 

2- In Chapter 17, we shall see that the answer is in general in the negative 
for the more general Riemannian spaces. 

3. The operation of putting a covariant index equal to a contravariant index in 
a tensor and summing over that common index is called contracticm. A contraction 
reduces the rank of a tensor by two: by one contravariant index and by one covari- 
ant index. 


Chapter 13 


1. The Navier-Stokes differential equations of hydrodynamics are discussed 
in Lamb’s Hydrodynamics. 


2. There are two viewpoints in hydrodynamics: one is the Eulerian point of 
view in terms of the Eulerian variables ; the other is the Lagrangean point of view 
in terms of the Lagrangean variables. For a non- viscous fluid, the Eulerian hydro- 
dynamical equations are the equations of motion of the fluid from the Eulerian point 
of view, and the Lagrangean hydrodynamical equations are the equations of motion 
Of the fluid from the Lagrangean point of view. The Eulerian hydrodynamical 
equations in rectangular coordinates y'^ are obtained by putting v = 0 in the Navier- 
Stokes equations 13-2. 


Eulerian Hydrodynamical Equations. 




1 Op 






0t/“ p by 

0 . 


If a*' are the rectangular coordinates of a fluid particle in the initial state of the 
fluid, and if the y^d'-, d?, a®, <) are the coordinates of the particle at time then the 
Lagrangean hydrodynamical equations are 

1 bp 


^iV bv^ 


p(2/S y\ y^) 


ba^ 


/ba^ " 

by^ by' 
bar da® 

bd^ da® 
d?/® dy® 
do‘- da® 


i=0. 


p ba^ 


= Po(a^, o®, o*). 


For a treatment of classical hydrodynamics, including treatments of the Eu- 
lerian differential equations and the Jjagrangean difTi'rcMitinl (‘ciiintions, the reader is 
referred to Lamb’s Hydrodynamics and to Webster’s Dynamics', (if. the references 
at the end of Part II. Ertel’s monograph on dynamic meteorology has some 
interesting remarks. 
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3. To obtain expansion 13 • 6 for the divergence we can proceed as follows. 
Since g, the determinant of the Euclidean metric tensor gry, is a relative scalar of 
weight two, it can be shown by the usual methods of obtaining covariant derivatives 
of absolute tensors that the co variant derivative g,i is given by 


Q,i 


ba;* 


- 2flrr 


cd 


But, in rectangular coordinates, g is unity and the Euclidean Christoffel symbols 
are zero. Hence g,i is zero in rectangular coordinates and consequently g,i « 0 in 
all coordinates. This means that 

blog-y/^ ' 


But the divergence of tt® is by definition 

bw® 


= 

ba^ 


+ 


so that by the above result we find the following equivalent expression for the divert 
gence 

„ 1 b(-\/gu^) 

= — 7 = • 

\/g 


If in particular u* = 


biA 

bx^ 


, where x®, x®) is a scalar field, we see that for 


this tt*' 

1 

“ -y/g e>x“ 

But in rectangular cartesian coordinates gr = 1, == 6®^, and so this scalar w,® 

reduces to the Laplacean in rectangular cartesian coordinates. Hence Laplnce's 
equation in general coordinates x* is given by 


1 

•^i/g bxi® 


0 when the unknown \^(x', x®, x®) is a scalar field. 


Chapter 14 

1. The fundamentals of a finite clastic deformation theory are not new. Kirch- 
hoff in 1852 made the firat systematic study, and E. and F. Cosserat in 1896 made 
an extensive investigation, of the subject. Ricci and Levi-Civita in 1900 made 
brief but important contributions to the applications of the tensor calculus to 
elasticity theory. L6on Drillouin in 1924 simplified and recast Cosserat’s treatment 
with the aid of the tensor calculus, and in 1937 F. I). Murnaghan, among several 
other authors, made contributions to the tensor theoretic ti’eatment of elasticity 
theory. 


Chapter 16 

1. One can consider strain differential invariants of order r, i.e., scalar fields 
that retain their forms as functions of the metric tensor gap, the strain tensor taPt 
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and the derivatives of cafi up to order r. Since successive covariant differentiation 
reduces to a corresponding order of partial differentiation in cartesian coordinates 
a strain differential invariant can be written exclusively in terms of tensor fields by 
merely replacing the derivatives of by corresponding covariant derivatives. The 
Michal-Thomas methods (see the 1927 papers of these authors) can be used to carry 
on some interesting researches on strain differential invariants. Strain different!^ 
tensors can also be considered. 

Examples of strain differential invariants of order one are 

B - and 9°^ ^ 

where ioep,y is the first covariant derivative of the strain tensor and where /i, 
/ 2 , and Iz are the three fundamental strain invariants. The question arises whether 
successive covariant differentiation of h, 1 2 , and combination with Qae would 
yield all the fundamental strain differential invariants. Clearly there exist no strain 
differential invariants for homogeneous strains. Note that the vanishing of R is the 
necessary and sufficient condition for<a homogeneous strain. 


Chapter 16 

1. The theory of complete systems of partial differential equations is discussed 
in Hedrick’s translation of Goursat’s Cours d’analyse, Vol. II, part II. 


Chapter 17 

1. For an account of Riemannian geometry, see Eisenhart’s Riemannian Geom- 
etry. This reference, of course, treats the (classical) finite dimensional Riemannian 
geometries. Infinite dimensional and dimensionless “Riemannian” geometries were 
first studied by A. D. Michal (see paper 2 under Michal in the references for Part II). 
The applications to vibrations of elastic media are now being studied (see papers 
4, 5, and 6 under Michal in the references). 

2. Several engineering applications of Lagrange’s equations of motion are to 
be found in the Kdrmdn and Biot book. 

3. We have seen that the Lagrangean equations of motion for a conservative 
dynamical system with no constraints and n degrees of freedom were 

(“) 

This means that the dynamical trajectories can be considered as curves in an n-dimen- 
( sional Riemannian space whose element of arc length ds is given by 

(6) ds® = gifiq^y • • - , 5")dg* d<fiy 

where the ga are the functions occurring in the kinetic energy T = hQijk^k^’ The 
differential equations of the dynamical curves are the second-order differential 
equations (o). 

These curves are not, in general, geodesics in the Riemannian space with arc 
lengths ds given by formula (6). The question then arises whether it is possible to 
define a Riemannian space whose geodesics are some of the curves whose differential 
equations are Lagrange’s equations of motion (a) for the given n.-degree dynamical 
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system. We shall show briefly that it is possible to define such a Riemannian 
space. To do this we must first show that the Lagrangean equations of motion have 
the energy integral 

r + F - C 

i.e., the sum of the hinetic and potential energies is a constant along any chosen 
namical trajectory. In fact, 

dL bLi ... bi . . 

Hence on using Lagrange’s equations of motion 

dL bL,.. d/dL\.. 

dt “ dt\,dg‘J^’ 

d/dL..\ 


We have then immediately T + V = C, a constant, since the kinetic potential 
L = T- V. 

Consider now the dynamical curves that correspond to any chosen energy constant 
C. We shall show that these particular dynamical curves are the geodesics of the 
n-dimensional Riemannian space whose element of arc length ds is given by 

(c) ds^ = 2(0- - V)gijiq\ dq^ d^\ 

where, as in (a), the gij are the coefficients in the kinetic energy T. For convenience in 
computation, let us define A = 2(C — F) and Oy = Agij so that (c) can be written 

(d) ds* = aij dq* d<ff. 

By definition 

,v Cofactor of an in a 

as 

a 

where a = the determinant of the aij. Hence, since factors out in the nu- 
merator and A" in the denominator respectively of we see that 


Let be the Christoffel symbols based on the metric tensor oij. Clearly 






^Ag„k bAg/g _ bAgjA; 

bq^ bg* bg^ ^ 


+ ds"”*/ 

By definition 

ds* - 2(C - V)gij dqi dq^ 

and hence along a d 3 niamical trajectory with energy constant C we have 


1 / bA 
lA V * bq^ 


. b A . bA " 


m- 


= C2(C - F)3* 


so that 
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along a dynamical trajectory with energy constant C. By elementary calculus we 
have therefore 


iff) 


ds dt * d^ 


d^q' 

dl^ 




dt dt 


The differential equations for the geodesics of the Riemannian space whose ds is 
given by (c) is 


Oi) 


ds® 3'^ds ds “ 


On employing results (e), (/), and (p) in (A), we get after obvious simplifications 


(O 

But 


di® dt dt 2A^ dt dt 


dq^ dc^ 


-A 


along a dynamical trajectory with energy constant C. Hence equations (i) reduce to 


0 ‘) 


di® ^^dl dt 


-9 


ptcr^ 




But these differential equations are another form of the Lagrangean differential 
equations of motion for our d3mamical system. We have thus proved that the 
geodesics of the RiemauTiian space with a ds given hy (c) are dynamical trajectories of 


the dynamical system with Tcin^tic energy T =» ^gij 


dt dt 


and potential energy V. 


By retracing our steps of proof, we can show that any dynamical trajectory with 
energy constant C, i.e., any curve that satisfies (j) with energy constant <7, can he con- 
sidered a geodesic in the Riemannian space with an element of arc length ds given in (c). 


Illustrative Example of a Shaft Carrying Four Disks. 

A shaft is fixed at one end and carries four disks at a distance I apart. If u is 
the moment of inertia of each disk and q^, q^, <fi the respective angular deflections 
of the four disks, then, if the shaft has a uniform torsional stiffness t, the kinetic 
and potential energies are given respectively by 

’■-C(s)’*(f)'*©' *(?)■] 

and 

y “ + it - tr ^it- tr ^it- ty-i- 


This is a conservative dynamical system of four degrees of freedom with no 
moving constraints. Hence the dynamical trajectories with total energy constant 
C can be represented as the geodesics of the four-dimensional Riemannian space 
whose element of arc length ds is given by 

ds8 » F{q\ gs, g^)C(dgi)* -h Wfi + + WY}. 

where 

yit, t‘s 3 ®, 2 *) “ + it- + it- ty + it- • 

Numerous other engineering examples can be given, some simpler and some 
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more sophisticated. For example, if the above shaft carries only tuoo disks, the 
dynamical trajectories with total energy constant C can be represented as the 
geodesics of the surface whose ds is given by 

- 2m jc - - 5‘)“3| C(<ia‘)= + 

A more sophisticated example is given by the symmetrical gyroscope; see the 
exercise at the end of Chapter 17. Here the dynamical trajectories with total 
energy constant C can be represented as the geodesics of the three-dimensional 
Riemannian space whose element of arc length ds is given by 

cZs® = 2i[C — Mgh cos + (-^ sin* ^ J cos* + 

2J cos ^ dq^ d<f + 


Chapter 18 


1. The conditions for an incompressible fluid and the continuity equations for 
a fluid flow state the invariance of two integrals: the integral for volume and the 
integral for mass, respectively. In other words, here we have two important ex- 
amples of integral invariants. For the theory of integral invariants and its modern 
generalizations, the reader is referred to paper 1 under Michal in the references. 
There the reader will find ample references to the earlier work on integral invariants 
by H. Poincar6, S. Lie, E. Cartan, and E. Goursat. 


2. We saw in Chapter 17 that the Riemann-Christoffel curvature tensor 
is a zero tensor in any Euclidean space and hence in our three-dimensional Eu- 
clidean space. Define the tensor (field) Riju by 

(tt) Rijkl — girBjJcU 

It is evident that 


(&) 


Rijkl * 0 


holds throughout our three-dimensional Euclidean space. It can be shown that 
there are only six independent equations in (6). Three of them are included in 

(c) Rocoffo •• Oj 

two of them are included in 

(d) RocPyo " Oj 

and the sixth one is given by 

(e) 1^1212 ™ 0. 

By straightforward calculation it can be shown that. 


I 


HctOpO 

I^otpyo 

R\2\2 


i: 1 

2 ~ 4^ da:® da:® 

PiSr;® — Toty;g 

*Ri2ii 4- (riir22 — FiaFai), 


where, it is to be recalled, 
surface a;" = C and = — 


a semicolon denotes covariant differentiation on the 
2 da::® ” 


The *Roefiy 5 stands for the curvature tensor on 


the surface. Define 

r 


bagCx\ a;*) 
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If we evaluate the last two sets of conditions ( /) on the surface = 0, we now see 
readily that the vanishing of the curvature tensor Rijki in the three-dimensional 
Euclidean space implies the following three sets of conditions : 

(ff\ ^ ^ 

2^ da^ 

(A) — ^ar/\P = 0 

(i) *RpctPy — bpfibocy “ bpyba^^ 

Equations (k) and (i) are the well-known Codazzi and Gauss equations of the 
surface = 0. (Cf. McConnell’s Applications of the Absolute Differential Calculus, 
p. 204, 1931.) 

Conversely it can be shown that equations (gr) in S-space and equations (Ji) and 
(i) over the surface a;° = 0 imply Rijki = 0 throughout the 3-space. But we shall not go 
into this matter any further. 

If we differentiate {g') with respect to afi and if in this result we eliminate the 
second derivatives by means of (^), we find 


(i) 


bi^Qac^ 1 bg^ ^gacp bpffy ^ 

2 dx® 2 

On differentiating the well-known identity 

Q\yg 


^ _ gU»g\y^£^ ^9acK ^ 


dx° da^ da;“ 


°x» 


we can solve for 


da^ 


and obtain 




dx" " " 5x° 

If we substitute this expression in (j), we evidently obtain 

£>®g«|8 


dx03 


0 . 


Hence the surface tensor components x®, x°) are quadratic expressions of x^- 

Indeed, if we write 


then 





» 


(k) gap =» Oceff + 2.boippfi + 

A glance at (g) shows that 

Cctfi “ a’^^^batrbpa- 

Hence, OotjS, bap, and Cap are respectively the tensor coefficients of the first, second, and 
third fundamental forms of the surface x° = 0. 

It can be shown that we may write (i) in the form 

CO ^^^Py^pot “ bppbay — bpybap, 

where 

€ap = Ot^VaP, ct *= 1 Oap I , ’flu =“ ^722 = 0, tji2 = 1» V2l = — 1» 

and 

K = l^^€P^*RpaPy = hU.a^^b,,f)^ - , 
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the total curvaiure (Gaussian curvature) of the surface a;® = 0. If we apply the tensor 
to (l)j we find 

(m) Cap - 2HhaS + “ 0 , 

where H is the mean curvature of the surface = 0, 

If Ri and Rz are the principal radii of curvature, it is well known in surface theory 
that 



With these relations, we can easily calculate 
For this purpose, we have to calculate the quantities 

if similar quantities involving Cfo, have been reduced with the help of (m). Now 

“ aa^^. 

Hence 

(n) if^rT^aayaps = 2a, « 2Ea. 

Applying to (Z)» we obtain 

(o) " 2K. 

If we substitute (n) and (o) in the determinant g, we find 

g.(l+2H^ + KWy-a{l+0(i+^J- 

Hence 

Thus, if we expand as a power series in a;'*, the series are convergent if [ a;® [ 
< Rm, where Rm is the minimum value of the principal radii of curvature of the 
surface = 0. The same is true of many other geometrical quantities derived 
from Qap and 
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ELEMENTARY APPLICATIONS 


of the Laplacean in curvilinear coordinates, the Navier-Stokes differ- 
ential equations take the following form in curvilinear coordinates x* 


(13-3) 


"ST 






b 

bx“ 



ba;“ 




bit 


+ (pw“).« = 0, 


where, as before, commas denote covariant differentiation based on the 
Euclidean Christoffel symbols x^, x^). If we expand the co- 

variant derivatives, we can write the Navier-Stokes differential equations 
in curvilinear coordinates as 


(13-4) 4 


bw* 


= ; 


H- 


/bn 




bw 


. bw'^ 

_j_ r 

~ ^ act 


du^ 


r>') - 1 


^ ia ^ /bw^ 


bw* 
bx 

^ + -X’S 

p OX 


.« d" rlaW 


Lm-) 


bp , b(pw“) 


bi 




+ = 0 . 


It is of interest in itself as well as in the above expansions to have 
equivalent expressions for the “ divergence in curvilinear coor- 
dinates. We have used the evident formula 


(13*5) 

but it can also he proved that^ 
(13-6) <« = 


l£\ =: 

•“ bx« 




1 ^(\/ g 

^/g bx“ 


where g is the determinant of the Euclidean metric tensor Formula 
13*6 is often more useful in the calculation of the divergence than 
formula 13*5. With the aid of formula 13.6 we see that the equation 
of continuity for a moving fluid takes the following foimi in curvilinear 
coordinoies x*: 


(13-7) 


bp(xS x\ x3, 0 b(pw“) , „^log g ^ 

d — T"" d- |pw“ - - „ = 0. 

ot bx“ bx“ 


Examples 

1. W^e saw in the last chapter that, if the x* are spherical polar 
coordinates, gw — 1, g^z = g^z — (x^)^ (sin x^)^, and all other 
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Qij *= 0. Hence 

and 


1 , 0 , 0 
0 , 0 
0 , 0 , {x^y(sin. x^y 


— (a:^)^(sin xP^y 


b log gr 4 




x^ 


^ logg 

brc^ 


= 2 cot a;^, 


^ loggr 


= 0 , 


The equation of continuity in spherical polar coordinates x^ then "becomes 


bi brc“ 


pi — r + cot x^ 


x^ 


) = o. 


2. As another exercise, we may take the x^ to be cylindrical polar 
coordinates so that 



Fig. 13- 1. 


ds^ = (dxy + ixy^dx^y + (dapy. 

Evidently 

g^ii = 1, 922 = (xy, gzz = 1, and all other gij = 0. 

Hence the determinant g = (a:')^ and we find readily 
the equation of continuity in cylindrical polar coordi- 
nates x^ 


^ b(p^“) ^ 

bi ba;“ ^x^ 


= 0 . 


Incidentally = 1, g^^ = 


{x}y 


g^^ = 1, and all other gr»? = 0, so that 


the Euclidean Christoffel symbols can easily be computed and found 

to be r 22 = — Tia = r|i = and all other Fk == 0. These ealcula- 

x^ 


tions for the Christoffel symbols are very much simpler than the cor- 
responding ones in Chapter 12 for spherical polar coordinates. 


Multiple-Point Tensor Fields, f 

The tensor fields that have been studied so far in this book have 
components that are functions of the coordinates of only one variable 
point in space. It is possible, however, to consider generalized tensor 
fields, called multiple-point tensor fields,! whose components depend 
on the coordinates of several points in space. Perhaps the simplest 
example of a two-point scalar field is the distance betweeri two points. 

t The first systematic research on multiple-point tensor fields was initiated by 
the writer many years ago. See A. D. Michal, Transactions of American Mathe- 
matical Society, vol. 29 (1927), pp. 612-646. 
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Let dQcif Xs) be the distance between two points having general 


coordinates (ic}, xf, xf) and 
(xl, Xz, xl) respectively. 
Each point may not neces- 
sarily be referred to the 
same coordinate system. 
For example, the xl may 






a; 2 = (a; J, «*,«») 


a;j= 

Fig. 13 - 2 . 

— j. ^ _ ^ 

be spherical coordinates while the X 2 may be cylindrical coordinates. 
Under transformation of coordinates 

(13-8) 4 = A(p^i, ^ 1 ) 

and 

(13-9) 4 


the components of distance transform by the rule 
(13-10) dCxi, X 2 ) = d(xi, X 2 ). 


Of course, if both points yi = (yl, yl, yl) and 2/2 = (yl, 2/i> 2/1) are referred 
to the same rectangular cartesian coordinate system, then the distance 
is given by the well-known formula 

(13 • 11) d(yt, 2/s) = j/g (i/i - y. 


Another simple example of a two-point tensor field is obtained as 
follows. Define 

(13-12) s(xit xz) = ICd(a;i, x^y, 

where d(xx, is the above two-point distance scalar. Obviously 
sixxf X 2 ) is also a two-point scalar field. Consider the partial derivatives 
of s(xx, X 2 ) with respect to the coordinates of the second point 


(13-13) 


Z>s(xi, xz) 


Now under transformations of coordinates 13-8 and 13-9 of the two 
points, we know that 

(13-14) s(5i, xz) = s{xi, Xz). 


Differentiating 13 • 14 we obtain 


(13-15) 


ds(:ci, Xz) 


Z)s(xi, Xz) ^X2 

bx2 ^xl 


This shows that the partial derivatives 13 • 13 are the components of 
a two-point tensor field: a covariant vector field with respect to the 
second point and a scalar field with respect to the first point. If both 
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points are referred to the same rectangular cartesian coordinaie system^ 
then from 13 • 11 we see that 


(13-16) 

If we define 
(13-17) 


2 / 2 ) 





^s(a:i, iCa) 


then obviously s^(a;i, x-f) is a two-point tensor field: a contravariant 
vector field with respect to the second point and a scalar field with 
respect to the first point. If both points are referred to the same 
rectangular coordinate system, then both two-point tensor fields 

^ and s*(a;i, x^ have the same components y\ — y\ in that 

0^2 

rectangular coordinate system. 

We shall have occasion to consider a two-point tensor field of rank 
two, a contravariant vector field with respect to each of the two points. 
As in the other examples, the two points need not be referred to the 
same coordinate system. The transformation law for the components 
of such a two-point tensor field are 

(13-18) i”(5i, ^ 2 ) = ^ 2 )^— 


A Two-Point Correlation Tensor Field in Turbulence. 

Let x^j x^, x^) be the contravariant velocity field of a fluid in 
motion. We shall denote the mean value of a function /(i) of the time 
t over the time interval (<o, ^i) by AfC/COJ* For example 

(13-19) r v(t, X\ di. 

t\ — toJto 

Clearly x^, x^, a:*)] is a contravariant vector field. Define a 

set of functions C^ixx, x^ of two points, Xx — (a:i, a:f, aj®) and x^ = 
(a4,a^, a|), by 

(13-20) 

. xxmt, x^n 

(Xi, X,, {g^(xdMZe<t, 

where ^®) is the Euclidean metric tensor in the general co- 
ordinates x\ It is evident that ^ 2 ) is a two-point tensor field of 

rank two, a contravariant vector field with respect to each of the two points; 
we shall call it the (two-point) correlation tensor field. If both points 
are referred to the same rectangular cartesian coordinate system, then 
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the correlation tensor field has components 


(13-21) 



in terms of the notations 

4 = 2/i). 



If we now assume that we are dealing with the special case of isotrojdc 
turhdme, the correlation tensor field 13-21 in rectangular coordinates 
simplifies still further and has components 

l M[u\u[\ 

3 M[{uY] 

on using the isotropic turbulence conditions that M[(u?)^] is inde- 
pendent of position and the index a, and equal, say, to 
Except for the numerical factor |, the above in rectangular coordi- 
nates is the correlation tensor used by Karman in isotropic turbulence. 
See his paper entitled “The Fundamentals of the Statistical Theory 
of Turbulence,” Journal of the AmmvMcal Sciences, vol. 4 (1937), 
pp. 131-138. 



CHAPTER 14 

APPLICATIONS OF THE TENSOR CALCULUS TO 

ELASTICITY THEORY 


Finite Deformation Theory of Elastic Media.^ 

One of the most natural and fruitful fields of application of the 
tensor calculus is to the deformation of media, elastic or otherwise. 
In the next three chapters we shall consider the fundamentals of the 
deformation of elastic media. We need not and shall not make the 
usual approximations of the classical (“infinitesimal’) theory in 
the general development of our subject. 

Consider a three-dimensional medium (a collection of point particles) 
in three-dimensional physical Euclidean space. We shall consider a 
deformation of the medium from the initial (unstrained) to its final 
(strained) position and obtain the strain tensor field under less stringent 
restrictions than those imposed in Chapter 10. 

Let (^a, ^a, ^a) be the curvilinear coordinates of a representative 
particle in an elastic medium, and let (ccS x^, x^) be the curvilinear 
coordinates of the representative particle after deformation. The 
deformation, a one-one point transformation A X, will be assumed 




given by differentiable functions 

(14-1) x^ = SK^a,‘^a,^a). 

It is convenient, and of some importance for the mathematical founda- 
tions, to assume that the representative unstrained particle A is repre- 
sented in a coordinate system not necessarily the same as the one in which 
the corresponding strained particle X is represented. For example, ^a, ^a, 
^a may be cylindrical coordinates while x^, x^, x^ are spherical polar 
coordinates. 

We shall adopt the following notational conventions with respect to 
one-point and two-point tensor fields. Tensor indices with respect to 
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transformation of coordinates of strained particles will be written to the 
righij and tensor indices with respect to transformation of coordinates 
of unstrained particles will be written to the left. For example, under a 
simultaneous transformation of the coordinates % of point A in the 
unstrained medium and of the coordinates of point X in the strained 
medium, 

(14-2) = - 

is a two-point tensor field of rank two, contravariant vector field with 
respect to point A and covariant vector field with respect to point X. 
In other words, under transforTnoiions of coordinates 

m.oN pa - ^<t>{^a, 2a, 

\x^ — X^, X^') 

in the unstrained and strained medium respectively, the two-point 
components *‘a,, undergo the transformation 

(14 • 4) ^a,g — a,p zr~ 

dx® c)“a 


Similarly 

(14-6) 

o*a 

is a two-point tensor field of rank two, covariant vector field with 
respect to point A and contravariant vector field with respect to 
point X. The relationships 

(14-6) (-o„) (.,*') - % = d: 

are clear, where 


rs - fir( = 0 if r 9^ s, 
) 1 

( = 1 if r = s. 


Let the initial and final squared elements of arc length in curvilinear 
coordinates be given respectively by 

(•iA.Y) ( dsQ — «/5^(<^)(d“a)(d^a), 

\ ds^ = dx“ dxP. 

Clearly, the initial and final squared elements in terms of the final 
coordinates x^ and initial coordinates *a respectively are 


(14-8) 

where 

(14-9) 


/ ds§ = h„^dx° dx^, 
\ds^ = pqk{d^a)(d^a)j 

(Kr = «/3C(“0(Vr), 


We are now in a position to write down the change produced by the 
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deformation in the squared arc element. In fact, in terms of the 
coordinates a;*, we have 

(14 • 10) ds^ — dsQ - 2e^^{x) dx** dxP, 

where 

(14*11) 

Similarly, in terms of the coordinates *a, we find 
(14 • 12) d^ - dsl » 2„^rjia) (d^a) (d^a), 

where 

(14*13) ccdnifl) = l(«/5^(«) - «/3c(a)). 

Clearly €„p(x) is a covariant tensor field of rank two in the strained” 
coordinates x* while is a covariant tensor field of rank two in the 

unstrained” coordinates *a. We shall call €^^( 3 ;) the Eulerian strain 
tensor and ccpvi^) "the Lagrangean strain tensor. The Eulerian strain 
tensor will often be referred to as the strain tensor. We have chosen this 
terminology in analogy with the two viewpoints in hydrodynamics 
represented respectively by the Eulerian and the Lagrangean differen- 
tial equations of motion. As an immediate consequence of formula 
14-10 we find the following fundamental result: A necessary and 
sufficient condition that the elastic deformation of the medium he a rigid 
motion (i.e., a degenerate deformation that merely displaces the medium 
in space with a preservation of distances between particles') is that the Eu- 
lerian strain tensor components be zero. Equivalently from 14- 12 we 
have: a necessary and sufficient condition for a rigid motion is the van- 
ishing of all the Lagrangean strain components. These results justify 
the use of the word “strain” in connection with the tensor fields 
and cts'n{o)- 


Strain Tensors in Rectangular Coordinates. 

If the same rectangular cartesian coordinate system is used for the 
description of both the initial and final positions of the elastic body, 
the Eulerian strain tensor reduces to 

(14-14) €^pix , a; , X ) - ^ — —y 


In terms of the usual notation (a, b, c) 
and (x, y, z) for the rectangular coordi- 
nates in the same coordinate system of 
the representative initial and final parti- 
cles respectively, we have 
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(14-15) 






€xs 


=2/1/ 


€j/z 




- ifi - - t—Y - VI 

■“ 21^ ” \dxj \dx/ ~ \dxJ J 

ir da da db db dc be"] 

2L dx dy dx dy dx b?/ J 
iF da da db db dc be"! 

2L dx dz dx dz dx bsJ 

4-(sr-(s)’-(s)’} 

iF ba da db db dc be"! 

2L dy dz dy dz dy dzj 

" 21^ “ \b^!/ ” \dz/ " VbJ J‘ 


■yxf 


zxy 


Similarly the Lagrangean strain tensor reduces to 

ba:^ ba:^ 


1/ ^ ba;^ ba:^ \ 

2\i^i 5^ 


(14-16) ^ £)“a d^o 

On denoting the “displacement vector” (really a two-point tensor 
field that was discussed in the previous chapter) {x — a, y — bf z — c) 
by (Uf V, w) we can rewrite 14 • 15 in the form 

du iF /bwY /bi/iVI 

ba; 2L\ba;/ \ba:/ ^ \ba:/ J 

du du 


€xx 




(14-17) I 


V— if- 

2\dv dx) ~ 2Lb: 

)4 


.dy 
l/du 


I 


2\b2 ba;, 


dv dv dw dw 

_j_ _l_ 

ba: dy dx dy dx dy 
du du dv dv dw bw”j 
ba: dz dx dz dx dz J 


Cjjxt 


= €■ 


■MXf 


tzyy 


tyz 


€zz ^ 


dv iF/b^Y fdv\^ /btuX'-^”] 
dy iLxbj// ^ \dy) ^ \dy ) 
l/dv bi«\ iF 

“2L 


dw 

dz 


If /duY 


I 




2L\b2 / 


+ 


du du dvdv dw biy"j 
dy dz dy dz dy dz J 




In the classical theory (the usual approximate theory) of elastic deforma- 
tions, the squares and products of the partial (lerivatives of u, v, and w are 
considered negligible. Hence, to the degree of approxhnation considered 
in the classical theory, 14 • 17 yields the following well-k^iown formulas 
for the strain tensor: 
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(14-18) 


£>w 


Ifbu 



“ Vx' 

— 


bx/ 

— ^yxy 

1/dw 




bz; 

~ 2\c)3 

. 

“ Caasj 

€yj/ 

by 

l/dv 

^w\ 



bw 

“ 2\5J 

" 

— €zyj 


bz 


Unstrained. 




Change in Volume under Elastic Deformation. 

Let us now return to our deformation theory without making the 
approximations of the classical theory. One of the first fundamental 
questions that arises is the manner in which volumes behave under a 

deformation of the medium. — 

The element of volume in 'N. ^ 

the unstrained medium is, / \ / 

by 12.4, / CD d 

(14 • 19) dUo = Vc d>a, d^a d^a, ^ 

where c = | q,^c | , the deter- Unstrained Strained 

minant of the «/ 3 C. Similarly, 14 • 3. 

(14-20) dV = y/g dx^ dx^ dx*, 

where g *= | g^^ | , the determinant of the g^^. 

In terms of the ‘‘strained variables’^ x*, we have 

(14-21) dVo = y/c I I dx^ dx^ dx®, 

I I - 

where “a «- is the determinant of the - — : - 


Strained 


Pig. 14-3. 


where “a «- is the determinant of the - — : - 
' ' bx* 

Now 

ds? = -K^cCob) 6^ a d'^a = /i«js(x) dx® dx^, 

where is given by formula 14 - 9. 

Evidently the determinant h of the is given by 

= C 1 “a.i 

From this 

(14-22) V* = Vc 1 “Oi 1 

and 

(14-23) dVo = y/h dx^ dx^ dx®. 

Formulas 14 - 20 and 14 • 23 imply that 


(14-24) 
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Define 

(14-25) 

and thus obtain 

(14-26) h^^gccXp 

by an application of the properties of g^p and see formulas 10-26. 
From 14 • 26 we compute 

(14-27) A = ? I h% \ , 

where I I is the determinant of the Kf. On using 14 ■ 24, we have 
immediately 

W “ ^TW\- 

To express this ratio in terms of the strain tensor we first recall the 
definition 14-11 and obtain 

Kp(x) = gccpipc) - 26«^(a^). 

On raising the indices with the aid of the g^, we evidently have 
(14-29) h%==5^-2e% 

where the mixed tensor field of rank two is defined by 
(14*30) e%ix) = g^-ixh^pix). 

On using 14-29 in 14 -28, we arrive at the important result that the 
ratio of the element of volume of a set of particles in the unstrained medium 
to the element of volume of the corresponding particles in the strained 
position is given in terms of the strain tensor eoep{.x) by means of the follow- 
iTig formula 

W ~ "^1 S|-26|(x) |. 

Notice that, with rigid motion of the medium, the strain tensor 
€ctp = 0 and hence — 0. Hence, from 14-31 and = 0, we see that 
a rigid motion preserves volumes. 

We have developed the fundamentals of elastic deformation and 
strain tensor in three-dimensional space. It is clear, however, that 
everything we have said can be taken over for the corresponding 
elastic deformations in the plane. In all the formulas, there will be 
two variables x^, x^, etc., the indices will have the range 1 to 2, and the 
consequent summations will go from 1 to 2. For example, formulas 
14-14, 14-15, 14-17, and 14-18 for the Eulerian strain tensor will be 
respectively as follows in elasticity theory in the plane. 
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(14-32) 


(14-33) 


€zx 




/ , ON 1/" 

, » ) = 

if ba ba bh bfeT ^ 
2L bx by bx bj/ J 

-i' -©■-©■} 


(14-34) 


Cjcr 


€xy 


ba; 2L Vba;/ \ba;/ J 

l/bw bv\ irbwbw bvdt;*! 
2\b2/ bx) A.bx by bx b^/J 




€yy 


bv 


5[©+(l)‘l 


where u = x ^ a and v - y - b, 

bu 


(14*35) 






bx 


l/bw bv 

“ rtl ^ + 


ijjy — 


2 \by 
bv 
by 


_ 

bx) ” 



Exercises 

1. Find the components of the Eulerian strain tensor eafiix^, a^, »“) when the 
deformation of the elastic body is a stretching whose equations are 

X* = A *a (4 is a constant greater than unity) 

where both x* and are rectangular coordinates referred to the same rectangular 
coordinate system. Discuss the change in volume elements. Work out the cor- 
responding problem in plane elasticity. 

2. Work out exercise 1 for a contraction so that the constant A is less than unity. 



CHAPTER 15 


HOMOGENEOUS AND ISOTROPIC STRAINS, STRAIN IN- 
VARIANTS, AND VARIATION OF STRAIN TENSOR 


Strain Invariants. 

If we expand the three-rowed determinant 
( 15 . 1 ) 
where 




we find 


b %- 2 e % I = 1 - 2/i + 4/2 - 8/3, 


h = sum of the principal two-rowed minors in the determinant 
(15. 2 ) A = I ej I , and I3 = A. 


A function f(jgii, gn, * • * , gss, €11, ci2, ■ * • , €33) of the Euclidean metric 
tensor and the strain tensor will he called a strain invariant^ if 
(a) it is a scalar field; (b) under all transformations of coordinates to 

( 15 - 3 ) f{gxi, ^12, • • •, gzzy in, ii2, * • *, €33) =/(s^ii, gvzj • • •, en, €12, 

* • €33), 

where the function /, on the left, is the same function of the and 
as it is, on the right, of the and e^p. 

We shall now prove that the three functions h, I2, and I3 occurring in 
the expansion of the determinant 15-1 are strain invariants. From the 
law of transformation of the mixed tensor field we readily get 

x^, 0) = x^, x^), 

from which follows that /i is a strain invariant on recalling that 
To prove that Iz is a strain invariant, we have by hypothesis 


i? (:^) = 6^(0:) 




On taking the determinant of corresponding sides, we obtain 




X 


bx^ 



€p 




5^® 

« 



But the product of the functional determinants is equal to unity. 
Hence | | = | \ , and I3 is a strain invariant. To prove that h 

is a strain invariant, we first observe that 



